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Abstract 

This paper surveys the application of geometric algebra to the physics 
of electrons. It first appeared in 1996 and is reproduced here with only 
minor modifications. Subjects covered include non-relativistic and rela- 
tivistic spinors, the Dirac equation, operators and monogenics, the Hydro- 
gen atom, propagators and scattering theory, spin precession, tunnelling 
times, spin measurement, multiparticle quantum mechanics, relativistic 
multiparticle wave equations, and semiclassical mechanics. 



Contents 

ll Introduction! 



E 



Sgacetime^^lgebraL 



2.1 The Spacetime SplTa 8 

^^^^pa^etim^^alculua 10 



l3 Spin ors an d the Dirac Equation! 12 



^^^^^^^gmor 



12 
17 



3.3 The Dirac Equation and Observable j 19 



!4 Ope rators^ Monogenics and the HydrogenAtom! 23 



4.1 Hamiltonian Form and the Non-F ^elativistic F^eductionl 23 



[cc 

4.2 Angular Ei genstates and Monogenic Functions! 27 

^"^Applkationa 30 



E-mail 



1 

^E-mail 
^E-mail 
*E-mail 



c.doran@mrao.cam.ac.uk 
a.n.lasenby@mrao.cam.ac.uk 
steve@mrao.cam.ac.uk 
a.d.challinor@mrao.cam.ac.uk 



Propagators and Scattering Theorvl 33 



5J_^^ro£a gajjx)iiaj^^ 34 

^^^^ 2inOT?otenMal^nT^roDagator^ . . 35 



5.3 Scattering Theorvl ... ... . . 35 



Plane Waves at Potential StepsI 38 



^J_^lgtchjii^_^^onditi^^ Waved 40 



^^^^atehin^m^^^ncsccn^ WavesI 42 

^^^^ p^^rocessi^^^^^^^^eil . . 44 



6^^^rujm£]^^ Wavea 46 

6.5 The Klein Parado^ 48 



TumiBlling_Tinies| 50 



7.1 Wavepacket Tunnellind . 50 



7.2 2-Diniensional Simulations! 55 



Sgh^^^easurement^ 55 



8.1 A Relativistic Model of a Spin Measurement! 57 

8'2"^^^vepa£kcr'simulatio^^ 59 



!9 The Multiparticle STA! 61 



9J^_2zPartjij£j^au]l.SteL^^ 64 



9.2 Com parison with the 'Causal' Approach toNon- Relativistic Spin Statea 69 
9'3"^£lativi^ti^''2'ParMd^^ 71 



^^^^lu]ti£aTti£lg_^^av e Eauationa 74 

9.5 The Pauh Principle! 76 



9.6 8-Dimcnsional Streamlines and Pauli Exclusion! 78 



EH 



^ 82 

10.1 Classical and Semiclassical Mechanic'3 82 



10.2 Grassmann Algebr 3 ........ 84 



!ll Conclusion^ 85 
!A The Spherical Monogenic Functions! 91 



2 



1 Introduction 



This paper surveys the apphcation of 'geometric algebra' to the physics of elec- 
trons. The mathematical ideas underlying geometric algebra were discovered 
jointly by Clifford T and Grassmann in the late 19th century. Their dis- 
coveries were made during a period in which mathematicians were uncovering 
many new algebraic structures (quaternions, matrices, groups, etc.) and the full 
potential of Clifford and Grassmann's work was lost as mathematicians concen- 
trated on its algebraic properties. This problem was exacerbated by Clifford's 
early death and Grassmann's lack of recognition during his lifetime. This paper 
is part of a concerted effort to repair the damage caused by this historical acci- 
dent. We firmly believe that geometric algebra is the simplest and most coherent 
language available for mathematical physics, and deserves to be understood and 
used by the physics and engineering communities. Geometric algebra provides a 
single, unified approach to a vast range of mathematical physics, and formulat- 
ing and solving a problem in geometric algebra invariably leeds to new physical 
insights. In the series of papers geometric algebra techniques were applied 

to number of areas of physics, including relativistic electrodynamics and Dirac 
theory. In this paper we extend aspects of that work to encompass a wider range 
of topics relevant to electron physics. We hope that the work presented here 
makes a convincing case for the use of geometric algebra in electron physics. 

The idea that Clifford algebra provides the framework for a unified language 
for physics has been advocated most strongly by Hestenes, who is largely re- 
sponsible for shaping the modern form of the subject. His contribution should 
be evident from the number and range of citations to his work that punctuate 
this paper. One of Hestenes' crucial insights is the role of geometric algebra 
in the design of mathematical physics 7 . Modern physicists are expected to 
command an understanding of a vast range of algebraic systems and techniques 
— a problem that gets progressively worst if one is interested in the theoretical 
side of the subject. A list of the some of the algebraic systems and techniques 
employed in modern theoretical physics (and especially particle physics) is given 
in Tabled Hestenes' point is that every one of the mathematical tools contained 
in Tabled can be expressed within geometric algebra, but the converse is not 
true. One would be hard-pressed to prove that two of the angles in an isosceles 
triangle are equal using spinor techniques, for example, but the proof is sim- 
ple in geometric algebra because it encompasses vector geometry. The work of 
physicists would be considerably simplified if, instead of being separately in- 
troduced to the techniques listed in Tabled they were first given a firm basis 
in geometric algebra. Then, when a new technique is needed, physicists can 
simply slot this into their existing knowledge of geometric algebra, rather than 
each new technique sitting on its own, unincorporated into a wider framework. 
This way, physicists are relieved of the burden of independently discovering the 
deeper organisational principle underlying our mathematics. Geometric algebra 
fulfills this task for them. 

In the course of this paper we will discuss a number of the algebraic systems 
listed in TableQ] and demonstrate precisely how they fit into the geometric alge- 
bra framework. However, the principle aim here is to discuss the application of 
geometric algebra to electron physics. These applications are limited essentially 
to physics in Minkowski spacetime, so we restrict our attention to the geometric 
algebra of spacetime — the spacetime algebra (HI. Our aim is twofold: to show 
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coordinate geometry 
complex analysis 



spinor calculus 
Grassmann algebra 
Berezin calculus 
differential forms 
twisters 

algebraic topology 



vector analysis 
tensor analysis 
Lie algebras 



Clifford algebra 



Table 1: 5*07716 algebraic systems employed in modern physics 



that spacetime algebra simplifies the study of the Dirac theory, and to show 
that the Dirac theory, once formulated in the spacetime algebra, is a powerful 
and flexible tool for the analysis of all aspects of electron physics — not just 
relativistic theory. Accordingly, this paper contains a mixture of formalism and 
applications. We begin with an introduction to the spacetime algebra (hence- 
forth the STA), concentrating on how the algebra of the STA is used to encode 
geometric ideas such as lines, planes and rotations. The introduction is designed 
to be self-contained for the purposes of this paper, and its length has been kept 
to a minimum. A list of references and further reading is given at the end of 
the introduction. 

In Sections 13 and 0] Pauli and Dirac column spinors, and the operators that 
act on them, arc formulated and manipulated in the STA. Once the STA for- 
mulation is achieved, matrices are eliminated from Dirac theory, and the Dirac 
equation can be studied and solved entirely within the real STA. A significant 
result of this work is that the unit imaginary of quantum mechanics is eliminated 
and replaced by a directed plane segment — a bivector. That it is possible to 
do this has many implications for the interpretation of quantum mechanics [HI . 

In SectionsEliniand[3we turn to issues related to the propagation, scattering 
and tunnelling of Dirac waves. Once the STA form is available, studying the 
properties of electrons via the Dirac theory is no more complicated than using 
the non-relativistic Pauli theory. Indeed, the first-order form of the Dirac theory 
makes some of the calculations easier than their non-relativistic counterparts. 
After establishing various results for the behaviour of Dirac waves at potential 
steps, we study the tunnelling of a wavepacket through a potential step. Indirect 
timing measurements for quantum-mechanical tunnelling are now available from 
photon experiments, so it is important to have a solid theoretical understanding 
of the process. We argue that standard quantum theory has so far failed to 
provide such an understanding, as it misses two crucial features underlying the 
tunnelling process. 

In Section |S| we give a relativistic treatment of a measurement made with a 
Stern-Gerlach apparatus on a fermion with zero charge and an anomalous mag- 
netic moment. As with tunnelling, it is shown that a disjoint set of outcomes 
is consistent with the causal evolution of a wavepacket implied by the Dirac 
equation. Wavepacket collapse is therefore not required to explain the results of 
experiment, as the uncertainty in the final result derives from the uncertainty 
present in the initial wavepacket. It is also argued that the standard quantum 
theory interpretation of the measurement performed by a Stern-Gerlach appa- 
ratus is unsatisfactory. In the STA, the anticommutation of the Pauli operators 
merely expresses the fact that they represent orthonormal vectors, so cannot 
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have any dynamical content. Accordingly, it should be possible to have simulta- 
neous knowledge of all three components of the spin vector, and we argue that a 
Stern-Gerlach apparatus is precisely what is needed to achieve this knowledge! 

Multiparticle quantum theory is considered in Section We introduce a 
new device for analysing multiparticle states — the multiparticle STA. This 
is constructed from a copy of the STA for each particle of interest. The re- 
sulting algebraic structure is enormously rich in its properties, and offers the 
possibility of a geometric understanding of relativistic multiparticle quantum 
physics. Some applications of the multiparticle STA are given here, including 
a relativistic treatment of the Pauli exclusion principle. The paper ends with a 
brief survey of some other applications of the STA to electron physics, followed 
by a summary of the main conclusions drawn from this paper. 

Summation convention and natural units {h = c = eq = 1) are employed 
throughout, except where explicitly stated. 

2 Spacetime Algebra 

'Spacetime algebra' is the name given to the geometric (Clifford) algebra gener- 
ated by Minkowski spacetime. In geometric algebra, vectors are equipped with 
a product that is associative and distributive over addition. This product has 
the distinguishing feature that the square of any vector in the algebra is a scalar. 
A simple re-arrangement of the expansion 

(a + bf = {a + b){a + b) ^ + [ab + ba) + b^ (2.1) 

yields 

ab + ba^{a + bf -a^ -b^, (2.2) 

from which it follows that the symmetrised product of any two vectors is also a 
scalar. We call this the inner product a -5, where 

a-b= \{ab + ba). (2.3) 

The remaining, antisymmetric part of the geometric product is called the outer 
product a A 6, where 

ahb= \{ab-ba). (2.4) 

The result of the outer product of two vectors is a bivector — a grade-2 object 
representing a segment of the plane swept out by the vectors a and b. 

On combining equations H2.3|l and (|2.4|l we see that the full geometric prod- 
uct of two vectors decomposes as 

ab = a-b + a/\b. (2.5) 

The essential feature of this product is that it mixes two different types of ob- 
ject: scalars and bivectors. One might now ask how the right-hand side of (|2.5|l 
is to be interpreted. The answer is that the addition implied by (|2.5() is that 
used when, for example, a real number is added to an imaginary number to form 
a complex number. We are all happy with the rules for manipulating complex 
numbers, and the rules for manipulating mixed-grade combinations are much 
the same But why should one be interested in the sum of a scalar and 
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a bivector? The reason is again the same as for complex numbers: algebraic 
manipulations arc simplified considerably by working with general mixed-grade 
elements (multivectors) instead of working independently with pure-grade ele- 
ments (scalars, vectors etc.). 

An example of how the geometric product of two vectors 1)2.5(1 is employed 
directly is in the description of rotations using geometric algebra. Suppose 
initially that the vector a is reflected in the hyperplane perpendicular to the 
unit vector n (n^ = !)• The result of this reflection is the vector 

a — 2a-nn = a — [an + na)n = —nan. (2-6) 

The form on the right-hand side is unique to geometric algebra, and is already 
an improvement on the usual formula on the left-hand side. If one now applies 
a second reflection in the hyperplane perpendicular to a second unit vector m 
the result is the vector 

—m{—nan)m — mna{mn)~. (2-7) 

The tilde on the right-hand side denotes the operation of reversion, which simply 
reverses the order of the vectors in any geometric product, 

{ab...cy=c...ba. (2.8) 

The combination of two reflections is a rotation in the plane specified by the 
two reflection axes. We therefore see that a rotation is performed by 

a i-> RaR (2.9) 

where 

R = mn. (2.10) 
The object R is called a rotor. It has the fundamental property that 

RR — mnnm = 1. (2-11) 

Equation ((2.9(1 provides a remarkably compact and efficient formulation for en- 
coding rotations. The formula for R ((2.10(1 shows that a rotor is formed from 
the geometric product of two unit vectors, so does indeed consist of the sum of a 
scalar and a bivector. A rotor can furthermore be written as the exponential of 
a bivector, R = ± exp{B/2), where the bivector encodes the plane in which the 
rotation is performed. This naturally generalises the complex representation of 
rotations frequently used in two dimensions. Rotors illustrate how mixed-grade 
objects are frequently employed as operators which act on other quantities in 
the algebra. The fact that both geometric objects and the operators that act on 
them are handled in a single unified framework is a central feature of geometric 
algebra. 

The above discussion applies to vector spaces of any dimension. We now 
turn to the case of specific interest, that of Minkowski spacetime. To make 
the discussion more concrete, we introduce a set of four basis vectors {7^}, 
= . . . 3, satisfying 

7m -71^ = VtJ^iy = diag(+ ). (2.12) 
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The vectors {7^} satisfy the same algebraic relations as Dirac's 7-niatrices, but 
they now form a set of four independent basis vectors for spacetime, not four 
components of a single vector in an internal 'spin space'. When manipulating 
(geometric) products of these vectors, one simply uses the rule that parallel 
vectors commute and orthogonal vectors anticommute. This result is clear im- 
mediately from equation H2.5|l . From the four vectors {7^^} we can construct a 
set of six basis elements for the space of bivectors: 

{7170, 7270, 7370, 7372, 7173, 7271 }■ (2-13) 

After the bivectors comes the space of grade-3 objects or trivectors. This space 
is four-dimensional and is spanned by the basis 

{737271: 707372, 7o7i 73, 7o727i}- (2-14) 

Finally, there is a single grade-4 element. This is called the pseudoscalar and is 
given the symbol i, so that 

i = 70717273. (2-15) 

The symbol i is used because the square of « is —1, but the pseudoscalar must not 
be confused with the unit scalar imaginary employed in quantum mechanics'^. 
The pseudoscalar z is a geometrically-significant entity and is responsible for the 
duality operation in the algebra. Furthermore, i aniicommutes with odd-grade 
elements (vectors and trivectors), and commutes only with even-grade elements. 
The full STA is spanned by the basis 

1> {7/x}, {o-fe, «crfc}, {ilp}, i, (2.16) 

where 

CTfc=7fc7o, fc= 1,2,3. (2.17) 

An arbitrary element of this algebra is called a multivector and, if desired, can 
be expanded in terms of the basis H2.16|l . Multivectors in which all elements have 
the same grade are usually written as Ar to show that A contains only grade-r 
components. Multivectors inherit an associative product from the geometric 
product of vectors, and the geometric product of a grade-r multivector Ar with 
a grade-s multivector Bs decomposes into 

ArBs - {AB)r+s + {AB)r+s-2 • ■ • + {AB)\r-s\- (2.18) 

The symbol {M)r denotes the projection onto the grade-r component of M . The 
projection onto the grade-0 (scalar) component of M is written as (M). The 
scalar part of a product of multivectors satisfies the cyclic reordering property 

{A...BC) = {CA...B). (2.19) 

The '•' and 'A' symbols are retained for the lowest-grade and highest-grade terms 
of the series (|2.18|l . so that 

Ar-Bs = {AB)\r-s\ (2.20) 
ArhBs - {AB)s+r, (2.21) 

^In more recent publications 1101 we have preferred the symbol / for the spacetime pseu- 
doscalar, and now recommend this ahead of the older i. 
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which are called the inner and outer (or exterior) products respectively. We also 
make use of the scalar product, defined by 

A^B={AB), (2.22) 

and the commutator product, defined by 

A-xB ^ \{AB - BA). (2.23) 

The associativity of the geometric product ensures that the commutator product 
satisfies the Jacobi identity 

Ax{By.C) + Bx{CxA) + Cy.{Ay.B) = 0. (2.24) 

When manipulating chains of products we employ the operator ordering con- 
vention that, in the absence of brackets, inner, outer and scalar products take 
precedence over geometric products. 

As an illustration of the working of these definitions, consider the inner 
product of a vector a with a bivector b he: 

a-{bAc) = {abhc)i 

— ^{abc — acb)i 

= a-bc — a-cb ~ ^(bac ~ cab)i. (2.25) 

The quantity bac — cab reverses to give minus itself, so cannot contain a vector 
part. We therefore obtain the result 

a-{bAc) ^ a-bc- a-cb, (2.26) 

which is useful in many applications. 

2.1 The Spacetime Split 

The three bivectors {ufc}, where Ufc — 7^70 (|2.17|l satisfy 

i(CTjCrfc + CTfcCTj) = -i(7i7fc + Iklj) = 5jk (2.27) 

and therefore generate the geometric algebra of three-dimensional Euclidean 
space 13 This is identified as the algebra for the rest-space relative to the 
timelike vector 70. The full algebra for this space is spanned by the set 

1, {^Tfc}, {zafc}, I, (2.28) 

which is identifiable as the even subalgebra of the full STA H2.1()(l . The iden- 
tification of the algebra of relative space with the even subalgebra of the STA 
simplifies the transition from relativistic quantities to observables in a given 
frame. It is apparent from 12.2711 that the algebra of the {ufc} is isomorphic 
to the algebra of the Pauli matrices. As with the {7;^}, the {ufe} are to be 
interpreted geometrically as spatial vectors (spacetime bivectors) and not as 
operators in an abstract spin space. It should be noted that the pseudoscalar 
employed in (|2.28|l is the same as that employed in spacetime, since 

0-20-3 = 7i7o727o737o = 7o7i7273 = «• (2.29) 
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The split of the six spacetime bivectors into relative vectors {cT/t} and relative 
bivectors {iak} is a frame-dependent operation — different observers determine 
different relative spaces. This fact is clearly illustrated using the Faraday bivec- 
tor F. The 'spacetime split' |H1E1 of F into the 70-system is made by separating 
F into parts which anticommute and commute with 79. Thus 

F = E + iB (2.30) 

where 

E^^iF^joFjo) (2.31) 

iB = ^{F + -foF^o)- (2.32) 

Both E and B are spatial vectors in the 70-frame, and iB is a spatial bivector. 
Equation (|2.3U|) decomposes F into separate electric and magnetic fields, and 
the explicit appearance of 70 in the formulae for E and B shows how this split 
is observer-dependent. Where required, relative (or spatial) vectors in the 70- 
system are written in bold type to record the fact that in the STA they are 
actually bivectors. This distinguishes them from spacetime vectors, which are 
left in normal type. No problems arise for the {cfc}, which are unambiguously 
spacetime bivectors, and so are left in normal type. 

When dealing with spatial problems it is useful to define an operation which 
distinguishes between spatial vectors (such as E) and spatial bivectors (such as 
iB). (Since both the {ak} and {iak} are spacetime bivectors, they behave the 
same under Lorentz-covariant operations.) The required operation is that of 
spatial reversion which, as it coincides with Hermitian conjugation for matrices, 
we denote with a dagger. We therefore define 

Aft ^ ^oAf7o, (2.33) 

so that, for example, 

F^ = E- iB. (2.34) 

The explicit appearance of 70 in the definition H2.33|) shows that spatial reversion 
is not a Lorentz-covariant operation. 

When working with purely spatial quantities, we often require that the dot 
and wedge operations drop down to their three-dimensional definitions. For ex- 
ample, given two spatial vectors a and 6, we would like aNb to denote the spatial 
bivector swept out by a and b. Accordingly we adopt the convention that, in 
expressions where both vectors are in bold type, the dot and wedge operations 
take their three-dimensional meaning. Whilst this convention may look clumsy, 
it is simple to use in practice and rarely causes any confusion. 

Spacetime vectors can also be decomposed by a spacetime split, this time 
resulting in a scalar and a relative vector. The spacetime split of the vector a 
is achieved via 

070 = a-7o + aA7o = ao + a, (2.35) 

so that flQ is a scalar (the 70-time component of a) and a is the relative spatial 
vector. For example, the 4-momentum p splits into 

p-fo = E + p (2.36) 
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where E is the energy in the 70 frame, and p is the S-momentuin. The definition 
of the relative vector (12.351) ensures that 



a-b= (a7o7o&) 

= {{ao + a){bo - b)} 

= aoba-a-b, (2.37) 
as required for the inner product in Minkowski spacetime. 

2.2 Spacetime Calculus 

The fundamental differential operator on spacetime is the derivative with respect 
to the position vector x. This is known as the vector derivative and is given the 
symbol V. The vector derivative is defined in terms of its directional derivatives, 
with the derivative in the a direction of a general multivector M defined by 

^,r/x M(x + ea)-M{x) 
a-VM(x) = lim — ^ 2.38) 

If we now introduce a set of four arbitrary basis vectors {e^}, with reciprocal 
vectors {e*^} defined by the equation ej-e'' — 6j, then the vector derivative 
assembles from the separate directional derivatives as 

V = e^ej-V. (2.39) 

This definition shows how V acts algebraically as a vector, as well as inheriting 
a calculus from its directional derivatives. 

As an explicit example, consider the {7;^} frame introduced above. In terms 
of this frame we can write the position vector x as x^^^, with x'^ = t, x^ = y. 
etc. and {x, y, z} a usual set of Cartesian components for the rest-frame of the 
7o vector. From the definition (|2.38|) it is clear that 



7.-V = ^ (2.40) 

which we abbreviate to 9^. From the definition (|2.39(l we can now write 

V = = 7°5t + 7'ax + 7'5y + fd, (2.41) 

which, in the standard matrix language of Dirac theory, is the operator that 
acts on Dirac spinors. It is not surprising, therefore, that the V operator should 
play a fundamental role in the STA formulation of the Dirac theory. What is 
less obvious is that the same operator should also play a fundamental role in the 
STA formulation of the Maxwell equations ^8 . In tensor notation, the Maxwell 
equations take the form 

a^F^'' = J^ = 0, (2.42) 

where [. . .] denotes total antisymmetrisation of the indices inside the bracket. 
On defining the bivector 

F^\F^''-f^A-1, (2.43) 
and the vector J = J^7p the equations p. 42(1 become 

V-F=J (2.44) 
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and 



VAF = 0. 



(2.45) 



But we can now utilise the geometric product to combine these separate equa- 
tions into the single equation 

VF = J (2.46) 

which contains all of the Maxwell equations. We see from (|2.46|) that the vector 
derivative plays a central role in Maxwell theory, as well as Dirac theory. The 
observation that the vector derivative is the sole differential operator required 
to formulate both Maxwell and Dirac theories is a fundamental insight afforded 
by the STA. Some consequences of this observation for propagator theory are 
discussed in [gj. 

The vector derivative acts on the object to its immediate right unless brack- 
ets are present, when it acts on everything in the brackets. Since the vector 
derivative does not commute with multivectors, it is useful to have a notation 
for when the derivative acts on a multivector to which it is not adjacent. We 
use overdots for this, so that in the expression VAB the V operator acts only 
on B. In terms of a frame of vectors we can write 

i/AB e^AeyVB. (2.47) 

The overdot notation provides a useful means for expressing Leibniz' rule via 

V{AB) = vAb + VAB. (2.48) 

The spacetime split of the vector derivative requires some care. We wish to 
retain the symbol V for the spatial vector derivative, so that 

V = akdk, fc = 1...3. (2.49) 

This definition of V is inconsistent with the definition H2.35|l , so for the vector 
derivative we have to remember that 



V70 = a* ~ V. (2.50) 

We conclude this introduction with some useful results concerning the vector 
derivative. We let the dimension of the space of interest be n, so that the results 
are applicable to both space and spacetime. The most basic results are that 

Va; = n (2.51) 

and that 

(2.52) 

where V' is an arbitrary multivector field. The latter result follows from the fact 
that partial derivatives commute. For a grade-r multivector Ar the following 
results are also useful: 

Wx-Ar^rAr (2.53) 
VxAAr ^ {n- r)Ar (2.54) 
VArX ^ {-iy{n-2r)Ar. (2.55) 
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More complicated results can be built up with the aid of Leibniz' rule, for 
example 

Vx'^ ^Vx-x + Vx-i = 2x. (2.56) 

This concludes our introduction to the spacetime algebra. Further details 
can be found in ^Space-Time Algebra' by Hestenes [S] and ^Clifford Algebra 
to Geometric Calculus' by Hestenes and Sobczyk |13j . The latter is a technical 
exposition of geometric algebra in general and does not deal directly with space- 
time physics. A number of papers contain useful introductory material including 
those by Hestenes ^lElEl and the series of papers OEEIEI written by three 
of the present authors. Baylis et al. and Void ^JE] have also written good 
introductory pieces, and the books 'New Foundations for Classical Mechanics' 
by Hestenes jll) and 'Multivectors and Clifford Algebras in Electrodynamics' by 
Jancewicz |18| provide useful background material. Further work can be found 
in the three conference proceedings |19ll2()ll?T] . though only a handful of papers 
are directly relevant to the work reviewed in this paper. Of greater interest are 
the proceedings of the conference entitled 'T/ie Electron' [221 1 which contains a 
number of papers dealing with the application of the STA to electron physics. 

3 Spinors and the Dirac Equation 

In this section we review how both the quantum states and matrix operators 
of the Pauli and Dirac theories can be formulated within the real STA. This 
approach to electron theory was initiated by Hestenes j^Sl 1^ and has grown 
steadily in popularity ever since. We start with a review of the single-electron 
Pauli theory and then proceed to the Dirac theory. Multiparticle states are 
considered in Sectional 

Before proceeding, it is necessary to explain what we mean by a spinor. The 
literature is replete with discussions about different types of spinors and their 
inter-relationships and transformation laws. This literature is highly mathe- 
matical, and is of very little relevance to electron physics. For our purposes, we 
define a spinor to be an element of a linear space which is closed under left-sided 
multiplication by a rotor. Thus spinors are acted on by rotor representations 
of the rotation group. With this in mind, we can proceed directly to study the 
spinors of relevance to physics. Further work relating to the material in this 
section is contained in 0]. 

3.1 Pauli Spinors 

We saw in Section [2.11 that the algebra of the Pauli matrices is precisely that 
of a set of three orthonormal vectors in space under the geometric product. So 
the Pauli matrices are simply a matrix representation of the geometric algebra 
of space. This observation opens up the possibility of eliminating matrices from 
the Pauli theory in favour of geometrically-significant quantities. But what of 
the operator action of the Pauli matrices on spinors? This too needs to be 
represented with the geometric algebra of space. To achieve this aim, we recall 
the standard representation for the Pauli matrices 
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The overhats distinguish these matrix operators from the {ak} vectors whose 
algebra they represent. The symbol i is reserved for the pseudoscalar, so the 
symbol j is used for the scalar unit imaginary employed in quantum theory. 
The {ak} operators act on 2-component complex spinors 



1*)^ (,;,;:]. (3.^) 

where ipi and '02 are complex numbers. Quantum states are written with bras 
and kets to distinguish them from STA multivectors. The set of \4>)^s form a 
two-dimensional complex vector space. To represent these states as multivectors 
in the STA we therefore need to find a four-dimensional (real) space on which 
the action of the {o^fc} operators can be replaced by operations involving the 
{tJfc} vectors. There are many ways to achieve this goal, but the simplest is 
to represent a spinor \ip) by an element of the even subalgebraof lj22H|)- This 
space is spanned by the set {1, iak} and the column spinor \ip) is placed in one- 
to-one correspondence with the (Pauli)-even multivector ip = 7o'07o through 
the identification 

m - [^y+fai) - = «° + (3.3) 



In particular, the spin-up and spin-down basis states become 

^ 1 (3.4) 



and 

^0 
1 



-ia2. (3.5) 



The action of the quantum operators {ak} and j is now replaced by the opera- 
tions 

CTfclV') ^ o-fe^cTa (/c = l,2,3) (3.6) 
and j\ip) ^ (3-7) 

Verifying these relations is a matter of routine computation; for example 

= I ■ 3 ^ -a^ + a^ia^ ~ a°ia2 + a^iai — aiipa3. (3.8) 

We have now achieved our aim. Every expression involving Pauli opera- 
tors and spinors has an equivalent form in the STA and all manipulations can 
be carried out using the properties of the {ak} vectors alone, with no need to 
introduce an explicit matrix representation. This is far more than just a theo- 
retical nicety. Not only is there considerable advantage in being able to perform 
the computations required in the Pauli theory without multiplying matrices to- 
gether, but abstract matrix algebraic manipulations arc replaced by relations of 
clear geometric significance. 
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Pauli Observables 

We now turn to a discussion of the observables associated with Pauh spinors. 
These show how the STA formulation requires a shift in our understanding of 
what constitutes scalar and vector observables at the quantum level. We first 
need to construct the STA form of the spinor inner product (V'l'/')- It is sufficient 
just to consider the real part of the inner product, which is given by 

^{m ^ ii^U), (3.9) 

so that, for example, 

= {a"f+a''a\ (3.10) 

(Note that no spatial integral is implied in our use of the bra-ket notation.) 
Since 

the full inner product becomes 

(V'l^) ^ (7A,</))s = (V'^0)--(V/</>Z(T3)*(73. (3.12) 

The right hand side projects out the {Ijicra} components from the geometric 
product tp^(j>. The result of this projection on a multivector A is written {A)s- 
For Pauli-even multivectors this projection has the simple form 

{A)s = U^-ifJsAias). (3.13) 

As an application of H3.12(l , consider the expectation value of the spin in the 
fc-direction, 

(V'lo-fcIV') ^ (ip^akipcra) ~ {'ip^ak'ipi)ia3. (3.14) 

Since tp^iakip reverses to give minus itself, it has zero scalar part. The right-hand 
side of (|3.14|l therefore reduces to 

(akipaaip^) ^ ak-{ipa3'ip^)v, (3.15) 

where (. . .)y denotes the relative vector component of the term in brackets. (This 
notation is required because (. . .)i would denote the spacetime vector part of 
the term in brackets.) The expression H3.15|l has a rather different interpretation 
in the STA to standard quantum mechanics — it is the crfc-component of the 
vector part of ipa^ip^ . As V'O'aV'^ is both Pauli-odd and Hermitian-symmetric it 
can contain only a relative vector part, so we define the spin- vector s by 

s^iPasip''. (3.16) 

(In fact, both spin and angular momentum are better viewed as spatial bivector 
quantities, so it is usually more convenient to work with is instead of s.) The 
STA approach thus enables us to work with a single vector s, whereas the 
operator/matrix theory treats only its individual components. We can apply 
a similar analysis to the momentum operator. The momentum density in the 
fc-dircction is given by 

{ipl-jdkW ^ -(V'Vfe-WifTa)- (V'Vfc-VV')itT3, (3.17) 
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in which the final term is a total divergence and so is ignored. Recombining 
with the {(7k} vectors, we find that the momentum vector field is given by 

p= -Vi^piasi}^). (3.18) 

It might appear that we have just played a harmless game by redefining 
various observables, but in fact something remarkable has happened. The spin- 
vector s and the momentum p are both legitimate {i.e. gauge-invariant) quan- 
tities constructed from the spinor ip. But standard quantum theory dictates 
that we cannot simultaneously measure all three components of s, whereas we 
can for p. The 'proof of this result is based on the non-commutativity of the 
{(T/c} operators. But, in the STA formulation, this lack of commutativity merely 
expresses the fact that the {ak} vectors are orthogonal — a fact of geometry, 
not of dynamics! Furthermore, given a spinor ^p there is certainly no difhculty 
in finding the vector s. So how then are we to interpret a spin measurement, 
as performed by a Stern-Gerlach apparatus for example? This problem will 
be treated in detail in Section |H1 but the conclusions are straightforward. A 
Stern-Gerlach apparatus is not a measuring device — it should really be viewed 
as a spin polariser. When a spinor wavepacket with arbitrary initial vector s 
enters a Stern-Gerlach apparatus, the wavepacket splits in two and the vector 
s rotates to align itself either parallel or anti-parallel to the B field. The two 
different alignments then separate into the two packets. Hence, in the final 
beams, the vector s has been polarised to point in a single direction. So, having 
passed through the apparatus, all three components of the spin- vector s are 
known - not just the component in the direction of the B field. This is a major 
conceptual shift, yet it is completely consistent with the standard predictions 
of quantum theory. Similar views have been expressed in the past by advocates 
of Bohm's 'causal' interpretation of quantum theory |2f)l 1271 1^ . However, the 
shift in interpretation described here is due solely to the new understanding of 
the role of the Pauli matrices which the STA affords. It does not require any 
of the additional ideas associated with Bohm's interpretation, such as quantum 
forces and quantum torques. 

Spinors and Rotations 

Further insights into the role of spinors in the Pauli theory are obtained by 
defining a scalar 

p = (3.19) 
so that the spinor ip can be decomposed into 

tP = p^/^R. (3.20) 

Here R is defined as 

R^p-^/^^ (3.21) 

and satisfies 

RR^ = 1. (3.22) 

In SectionElwe saw that rotors, which act double-sidedly to generate rotations, 
satisfy equation (|3.22|) . It is not hard to show that, in three dimensions, all even 
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s = ■(/'O'aV'^ 



Table 2: Summary of the main results for the STA representation of Pauli 
spinors 



quantities satisfying H^122|l are rotors. It follows from (|3.2Q(I that the spin- vector 
s can now be written as 

s^pRa^R^ (3.23) 

which demonstrates that the double-sided construction of the expectation value 
in equation H3.15|l contains an instruction to rotate the fixed 173 axis into the 
spin direction and dilate it. The decomposition of the spinor ■0 into a density 
term p and a rotor R suggests that a deeper substructure underlies the Pauli 
theory. This is a subject which has been frequently discussed by Hestenes |^ 
1291 130L ini . As an example of the insights afforded by this decomposition, one 
can now 'explain' why spinors transform single-sidedly under rotations. If the 
vector s is to be rotated to a new vector RqsRq then, according to the rotor 
group combination law, R must transform to RqR. This induces the spinor 
transformation law 

^ i?oV' (3.24) 
which is the STA equivalent of the quantum transformation law 

j 

\tp) ^ exp{-enkak}\^j) (3.25) 

where {uk} are the components of a unit vector. 

We can also now see why the presence of the 173 vector on the right-hand side 
of the spinor ip does not break rotational invariance. All rotations are performed 
by left-multiplication by a rotor, so the spinor effectively shields the CT3 on 
the right from the transformation. There is a strong analogy with rigid-body 
mechanics in this observation, which has been discussed by Hestenes ^ |^ . 
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Similar ideas have also been pursued by Dewdney, Holland and Kyprianidis |33 
I33| . We shall see in the next section that this analogy extends to the Dirac 
theory. The main results of this section are summarised in Table |2] 



3.2 Dirac Spinors 

The procedures developed for Pauli spinors extend simply to Dirac spinors. 
Again, we seek to represent complex column spinors, and the matrix operators 
acting on them, by multivectors and functions in the STA. Dirac spinors are four- 
component complex entities, so must be represented by objects containing 8 real 
degrees of freedom. The representation that turns out to be most convenient 
for applications is via the 8-dimensional even subalgebra of the STA |^ . If 
one recalls from Section ITTl that the even subalgebra of the STA is isomorphic 
to the Pauli algebra, we see that what is required is a map between column 
spinors and elements of the Pauli algebra. To construct such a map we begin 
with the 7-matrices in the standard Dirac-Pauli representation j35| , 

70= (n °r)' I'fe^li'. n) ^nd 75 = ( ^ ^ ) , (3.26) 



-/y ' \ak J \I 

where 75 = 7^ = — i7o7i7273 and / is the 2x2 identity matrix. A Dirac 
column spinor \^) is placed in one-to-one correspondence with an 8-component 
even element of the STA via P^l^ 



—a} + ja^ 



V' = a° + a''icrfc + (6" + 6'=icrfe)cr3. (3.27) 



With the spinor ji/)) now replaced by an even multivector, the action of the 
operators {7^^, 75,^} becomes 

^ 7^^-70 (m = 0,...,3) (3.28) 

^ ^10^ (3.29) 

75 IV') ^ V'o-a- (3.30) 

To verify these relations, we note that the map (|3.27|l can be written more 
concisely as 



IV-) = (^'iJj ^ V = '/' + W3, (3.31) 

where |(^) and I77) are two-component spinors, and </) and r\ are their Pauli-even 
equivalents, as defined by the map H3.3|l . We can now see, for example, that 

IkW) = oljl^^) ^ ^f^fcWs + crfe0 = 7/c(0 + W3)7o, (3.32) 

as required. The map (|3.31() shows that the split between the 'large' and 'small' 
components of the column spinor ji/)) is equivalent to splitting -0 into Pauli-even 
and Pauli-odd terms in the STA. 
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Alternative Representations 

All algebraic manipulations can be performed in the STA without ever intro- 
ducing a matrix representation, so equations H3.27|l and l|3.30|l achieve a map 
to a representation-free language. However, the explicit map H3.27(l between 
the components of a Dirac spinor and the multivector ijj is only relevant to the 
Dirac-Pauli matrix representation. A different matrix representation requires 
a different map so that that the effect of the matrix operators is still given 
by H3.30|l . The relevant map is easy to construct given the unitary matrix S 
which transforms between the matrix representations via 

% = S%S-\ (3.33) 

The corresponding spinor transformation is {ip) ^— > Slip), and the map is con- 
structed by transforming the column spinor j^)' in the new representation back 
to a Dirac-Pauli spinor S^^p)' . The spinor S^ip)' is then mapped into the 
STA in the usual way H3.27|l . As an example, consider the Weyl representation 
defined by the matrices 

,i ^ ( = {I i;) ■ (»^) 

The Weyl representation is obtained from the Dirac-Pauli representation by the 
unitary matrix 



A spinor in the Weyl representation is written as 



I0)'=(KV (3.36) 



where \x) and |^) are 2-component spinors. Acting on with gives 

^^i^)' = :).fK:!SV (3.37) 



V2 y\x) + \fj) 

Using equation (|3.27|) . this spinor is mapped onto the even element 

^^1^)' - 71 (1^) ; II) ^ = ^75 (1 + '^3) - ^^(1 - aa), (3.38) 

where x a^id fj are the Pauli-even equivalents of the 2-component complex 
spinors \x) and \fj), as defined by equation 1)3.3(1 . The even multivector 

V^ = ^72(l + ^3)-r7^(l-f73) (3.39) 

is therefore our STA version of the column spinor 

W-Oii], (3.40) 



where \ip)' is acted on by matrices in the Weyl representation. As a check, we 
observe that 

loW = (Z|jjj) -^7l(l + ^3) + X7f (1 - '^3) = 7oV'7o (3.41) 
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and 

TfelV^) = ^"^l"!^^^ ^ -C^fcW3^(l + cr3)-crfeXCr3^(l-cr3) = TfcV'To- (3.42) 

(Here we have used equation H3.7() and the fact that 70 commutes with all Pauli- 
even elements.) The map (I3.38|) does indeed have the required properties. 

While our procedure ensures that the action of the {7^, 75} matrix operators 
is always given by l|3.3()|l . the same is not true of the operation of complex 
conjugation. Complex conjugation is a representation-dependent operation, so 
the STA versions can be different for different representations. For example, 
complex conjugation in the Dirac-Pauli and Weyl representations is given by 

It/))* ^ -72V'72, (3.43) 

whereas in the Majorana representation complex conjugation leads to the STA 
operation ^ 

l^)Maj ^^^2. (3.44) 

Rather than think of H3.43|l and H3.44(l as different representations of the same 
operation, however, it is simpler to view them as distinct STA operations that 
can be performed on the multivector ^. 

3.3 The Dirac Equation and Observables 

As a simple application of H3.27() and (|3.30() . consider the Dirac equation 

rUd.-eAM^Mi^)- (3-45) 

The STA version of this equation is, after postmultiplication by 70, 

Vipia^ — eAijj — rmpjo, (3.46) 

where V = 7'^ 9^ is the spacetime vector derivative (j2.41|l . The STA form of 
the Dirac equation (|3.46|l was first discovered by Hestenes , and has been dis- 
cussed by many authors since; see, for example, references (36., ,38,. 39, 4Q. 4fj. 
The translation scheme described here is direct and unambiguous and the result- 
ing equation is both coordinate-free and representation-free. In manipulating 
equation (|3.46(l one needs only the algebraic rules for multiplying spacetime 
multivectors, and the equation can be solved completely without ever introduc- 
ing a matrix representation. Stripped of the dependence on a matrix represen- 
tation, equation H3.46|l expresses the intrinsic geometric content of the Dirac 
equation. 

In order to discuss the observables of the Dirac theory, we must first consider 
the spinor inner product. It is necessary at this point to distinguish between 
the Hermitian and Dirac adjoint. These are written as 

(V^l - Dirac adjoint , . 

— Hermitian adjoint, 

which are represented in the STA as follows, 
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One can see clearly from these definitions tliat the Dirac adjoint is Lorentz- 
invariant, whereas the Hcrmitian adjoint requires singling out a preferred time- 
like vector. 

The inner product is handled as in equation (|3.12(l . so that 

{^\4>) ^ {U)-{U^<^:i)^^z = {U)s, (3.49) 

which is also easily verified by direct calculation. By utilising H3.49|l the STA 
forms of the Dirac spinor bilinear covariants |37| are readily found. For example, 

(V'ITmIV') ^ (V'7mV'7o) - ('07MV'«73)jcr3 = 7Ar('07oV')i (3-50) 

identifies the 'observable' as the 7^-component of the vector {ipjo'ip)i. Since the 
quantity tpjotp is odd and reverse-symmetric it can only contain a vector part, 
so we can define the frame-free vector J by 

J = ^70^- (3.51) 

The spinor ip has a Lorentz-invariant decomposition which generalises the 
decomposition of Pauli spinors into a rotation and a density factor (|3.20(l . Since 
iptp is even and reverses to give itself, it contains only scalar and pseudoscalar 
terms. We can therefore define 

pe"^ = ^4, (3.52) 
where both p and /3 are scalars. Assuming that p ^ 0, ip can now be written as 

^ = p^^e'P/^R (3.53) 

where 

R ^ (pe''')-i/V- (3.54) 

The even multivector R satisfies RR = 1 and therefore defines a spacetime rotor. 
The current J HH.51|I can now be written as 

J ^ pv (3.55) 

where 

V = RjqR. (3.56) 

The remaining bilinear covariants can be analysed likewise, and the results are 
summarised in Table |3| The final column of this Table employs the quantities 

s = ipl3'4'j and S — ipia^ip. (3.57) 

Double-sided application of i? on a vector a produces a Lorentz transfor- 
mation 3 . The full Dirac spinor -0 therefore contains an instruction to rotate 
the fixed {7^} frame into the frame of observables. The analogy with rigid- 
body dynamics first encountered in Section 13.11 with Pauli spinors therefore 
extends to the relativistic theory. In particular, the unit vector v (I3.56|l is both 
future-pointing and timelike and has been interpreted as defining an electron 
velocity j^O] (see also the critical discussion in [SI). The factor' appearing 
in the decomposition of ip (|3.52|) has also been the subject of much discus- 
sion |S1 lUl I39L I4U| since, for free-particle states, /3 determines the ratio of particle 
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Bivector 
Pseudovcctor 
Pseudoscalar 


(V^IV) 

(V']7m75|V') 

(V'|j75|^) 


7m- ("^70-0) 
(7m^7i^)-(V'«_o-3V') 
7m-('/'73V') 

{ip-ipi) 


p cos/3 

V'7oV^= J 
ipia^ip — S 

ipla'ip = s 
—p sin/3 



Table 3: Bilinear covariants in the Dirac theory. 



to anti-particle solutions. It remains unclear whether this idea extends usefully 
to interacting systems. 

In Section ini we argued that, for Pauli spinors, any dynamical consequences 
derived from the algebraic properties of the Pauli matrices were questionable, 
since the algebra of the Pauli matrices merely expresses the geometrical rela- 
tions between a set of orthonormal vectors in space. Precisely the same is true 
of any consequences inferred from the properties of the Dirac matrices. This 
observation has the happy consequence of removing one particularly prevalent 
piece of nonsense — that the observed velocity of an electron must be the speed 
of hght The 'proof of this result is based on the idea that the velocity 

operator in the fc-direction is the 'jklo matrix. Since the square of this matrix 
is 1, its eigenvalues must be ±1. But in the STA the fact that the square of the 
{7^} matrices is ±1 merely expresses the fact that they form an orthonormal 
basis. This cannot possibly have any observational consequences. To the extent 
that one can talk about a velocity in the Dirac theory, the relevant observable 
must be defined in terms of the current J. The {7^^} vectors play no other role 
than to pick out the components of this current in a particular frame. The 
shift from viewing the {7^^} as operators to viewing them as an arbitrary, fixed 
frame is seen clearly in the definition of the current J (|3.51|l . In this expres- 
sion it is now the ip that 'operates' to align the 70 vector with the observable 
current. Since ip transforms single-sidedly under rotations, the fixed initial 70- 
vector is never affected by the rotor and its presence does not violate Lorentz 
invariance 0]. 

We end this subsection by briefly listing how the C, P and T symmetries 
are handled in the STA. Following the conventions of Bjorken & Drell we 
find that 

P\ip) ^ 70-0(^)70 

(710) ^ 0(71 (3.58) 

f|0;) ^ i7o0(-i)7i, 

where x = "fox^o is (minus) a reflection of x in the timelike 70 axis. The 
combined CPT symmetry corresponds to 

ip ^ -iip{-x) (3.59) 

so that CPT symmetry does not require singling out a preferred timelike vector. 
A more complete discussion of the symmetries and conserved quantities of the 
Dirac theory from the STA viewpoint is given in [5]. There the 'multivector 
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derivative' was advocated as a valuable tool for extracting conserved quantities 
from Lagrangians. 



Plane-Wave States 

In most applications of the Dirac theory, the external fields applied to the elec- 
tron define a rest-frame, which is taken to be the 70-frame. The rotor R then 
decomposes relative to the 70 vector into a boost L and a rotation $, 

R = i$, (3.60) 

where 

= L (3.61) 

$^ = $ (3.62) 

and LL = $$ = 1. A positive-energy plane- wave state is defined by 

^ = ^oe"^''='^'^ (3.63) 

where ipo is a constant spinor. From the Dirac equation 13.46|l with A = 0, it 
follows that tpo satisfies 

pV'o = rmpojo- (3.64) 
Postmultiplying by "00 we see that 

piptp = mJ (3.65) 

from which it follows that exp(i/3) = ±1. Since p has positive energy we must 
take the positive solution {(3 — 0). It follows that is just a rotor with a 
normalisation constant. The boost L determines the momentum by 

p — mLjoL = ?7iL^7o, (3.66) 

which is solved by 

r / r~ E + m + p 

L = ^Pl^lm = (3.67) 
yj2m{hj + m) 

where 

P7o==-B+p. (3.68) 

The Pauli rotor <i> determines the 'comoving' spin bivector ^iiTa^. This is 
boosted by L to give the spin S as seen in the laboratory frame. A $173$ gives 
the relative spin in the rest-frame of the particle, we refer to this as the 'rest- 
spin'. In Section f6.3l we show that the rest-spin is equivalent to the 'polarisation' 
vector defined in the traditional matrix formulation. 

Negative energy solutions are constructed in a similar manner, but with an 
additional factor of i or CT3 on the right (the choice of which to use is simply a 
choice of phase) . The usual positive- and negative-energy basis states employed 
in scattering theory are (following the conventions of Itzykson & Zuber |37l 
Section 2-2]) 

positive energy V^^^a^) = Wr(p)e"""'P '' (3.69) 
negative energy V^^^a;) = u^(p)e""=^P-'' (3.70) 
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with 



Ur{p) = L{p)Xr 
Vr{p) = L{p)XrCr3- 



(3.71) 
(3.72) 



Here L{p) is given by equation (|3.t)7|l and Xr — {1, are spin basis states. 

The decomposition into a boost and a rotor turns out to be very useful in 
scattering theory, as is demonstrated in Sectional 

The main results for Dirac operators and spinors are summarised in Table 01 

4 Operators, Monogenics and the Hydrogen 
Atom 

So far, we have seen how the STA enables us to formulate the Dirac equation 
entirely in the real geometric algebra of spacetime. In so doing, one might worry 
that contact has been lost with traditional, operator-based techniques, but in 
fact this is not the case. Operator techniques are easily handled within the STA, 
and the use of a coordinate-free language greatly simplifies manipulations. The 
STA furthermore provides a sharper distinction between the roles of scalar and 
vector operators. 

This section begins by constructing a Hamiltonian form of the Dirac equa- 
tion. The standard split into even and odd operators then enables a smooth 
transition to the non-relativistic Pauli theory. We next study central fields and 
construct angular-momentum operators that commute with the Hamiltonian. 
These lead naturally to the construction of the spherical monogenics, a basis 
set of orthogonal eigenfunctions of the angular-momentum operators. We finally 
apply these techniques to two problems — the Hydrogen atom, and the Dirac 
'oscillator'. 

4.1 Hamiltonian Form and the Non-Relativistic Reduc- 



The problem of how to best formulate operator techniques within the STA is 
really little more than a question of finding a good notation. We could of course 
borrow the traditional Dirac 'bra-kef notation, but we have already seen that 
the bilinear covariants are better handled without it. It is easier instead to just 
juxtapose the operator and the wavefunction on which it acts. But we saw in 
Section O that the STA operators often act double-sidedly on the spinor ip. This 
is not a problem, as the only permitted right-sided operations are multiplication 
by 7o or itT3, and these operations commute. Our notation can therefore safely 
suppress these right-sided multiplications and lump all operations on the left. 
The overhat notation is useful to achieve this and we define 



It should be borne in mind that all operations are now defined in the STA, so 
the 7^ are not intended to be matrix operators, as they were in Section lO 

It is also useful to have a symbol for the operation of right-sided multiplica- 
tion by ia^. The symbol j carries the correct connotations of an operator that 



tion 



7m V' = 7m "070- 



(4.1) 
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Equation 



VV'io'3 — eAip = mipjo 



Observables 



S = tpia^'ip s = f/'TsV' 



Plane- Wave 
States 



L{p) = (p7o + m)/y/2m{E + m) 



Table 4: Summary of the main results for the STA representation of Dirac 
spinors. The matrices and spinor equivalence are for the Dirac-Pauli represen- 
tation. The spinor equivalences for other representations are constructed via 
the method outlined in the text. 
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commutes with all others and squares to —1, and we define 



jqp = il^ia^. (4-2) 
The Dirac equation H3.46|l can now be written in the 'operator' form 

jVip — eAip — rmp. (4-3) 

where 

V')/; = VtAto, and A-tp = Aip^fQ. (4.4) 

Writing the Dirac equation in the form (|4.3|l docs not add anything new, but 
does confirm that we have an efficient notation for handling operators in the 
STA. 

In many applications we require a Hamiltonian form of the Dirac equation. 
To express the Dirac equation ()3.46(l in Hamiltonian form we simply multiply 
from the left by 70. The resulting equation, with the dimensional constants 
temporarily put back in, is 

jhdtip ~ cpip + eVtp — ceAip + mc^ip (4-5) 

where 

p^/; = -jhVip (4.6) 

■0 = Jo^lo (4.7) 

and 7oA -cA. (4.8) 

Choosing a Hamiltonian is a non-covariant operation, since it picks out a pre- 
ferred timelike direction. The Hamiltonian relative to the 70 direction is the 
operator on the right-hand side of equation 1)4. 5|l . We write this operator with 
the symbol H. 



The Pauli Equation 

As a first application, we consider the non-relativistic reduction of the Dirac 
equation. In most modern texts, the non-relativistic approximation is carried 
out via the Foldy- Wouthuysen transformation [SH EZ) • Whilst the theoretical 
motivation for this transformation is clear, it has the defect that the wavefunc- 
tion is transformed by a unitary operator which is very hard to calculate in all 
but the simplest cases. A simpler approach, dating back to Feynman 021 j is 
to separate out the fast-oscillating component of the waves and then split into 
separate equations for the Pauli-cven and Pauli-odd components of ip. Thus we 
write (with h — I and the factors of c kept in) 

V; = (0 + 77)e-^'"3"^='* (4.9) 

where (j) = <p and fj — —77. The Dirac equation (|4.5(l now splits into the two 
equations 

£(l)~cOr]^0 (4.10) 
{£ + 2mc^)T]- cO<j) = 0, (4.11) 
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where 



£0 ^ [jdt - eV)<t> (4.12) 
0(/) = (p - eA)0. (4.13) 

The formal solution to the second equation H4.11|l is 

^ = ^{\ + :r^) 0<i>, (4.14) 

where the inverse on the right-hand side is understood to denote a power series. 
The power series is well-defined in the non-relativistic limit as the £ operator is 
of the order of the non-relativistic energy. The remaining equation for is 

^^-|^(l-^ + ->^ = 0> (4.15) 

which can be expanded out to the desired order of magnitude. There is little 
point in going beyond the first relativistic correction, so we approximate (|4.15() 

by 

4m^c^ 2m 

We seek an equation of the form £0 — Hcj), where H is the non-relativistic 
Hamiltonian. We therefore need to replace the 0£0 term in equation (|4.1t)|) by 
a term that does not involve £. To do so we would like to utilise the approximate 
result that 

S^f^^c^, (4.17) 

but we cannot use this result directly in the 0£0 term since the £ does not 
operate directly on 0. Instead we employ the operator rearrangement 

20£0 ^[0,[£,0]]+£0^ + 0^£ (4.18) 

to write equation H4.16|l in the form 

£0^ + 0^£ 1 



We can now make use of H4.17II to write 



£0^<j)KO'^£(j)K ^ + 0(c-^) (4.20) 
2m 

and so approximate H4.16|l by 

£cj,^^<l^--^[0,[£,0]]q^~-^q^, (4.21) 
which is valid to order c~^. The commutators are easily evaluated, for example 
[£, O] = ~je{dt A + VV) = jeE. (4.22) 
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There are no time derivatives left in this commutator, so we do achieve a sensible 
non-relativistic Hamiltonian. The full commutator required in equation (|4.21() 
is 

[0,[£,0]]^[~3V -eA.jeE] 

^ {eVE)~2eEAV -2je^AAE (4.23) 

in which the STA formulation ensures that we are manipulating spatial vectors, 
rather than performing abstract matrix manipulations. 

The various operators H4.12|l . (|4.13|) and l|4.23|l can now be fed into equa- 
tion H4.21|l to yield the STA form of the Pauli equation 

- tAt^H'^E - 2EAVU - 2e^AAE(j>ia3], (4.24) 

which is valid to 0(c~^). (We have assumed that \A\ ~ c^^ to replace the O'* 
term by p^.) Using the translation scheme of Table El it is straightforward to 
check that equation (|4.24l) is the same as that found in standard texts In 
the standard approach, the geometric product in the 'VE term (|4.24|l is split 
into a 'spin-orbit' term VA-E and the 'Darwin' term 'V E. The STA approach 
reveals that these terms arise from a single source. 

A similar approximation scheme can be adopted for the observables of the 
Dirac theory. For example the current, V'ToV'j has a three-vector part 

J = (V'7oV')A7o = + ri4>\ (4.25) 
which is approximated to first-order by 

J« -l((V0z(730^)„~ (4.26) 

TO 

Not all applications of the Pauli theory correctly identify H4.26|l as the conserved 
current in the Pauli theory — an inconsistency first noted by Hestenes and 
Gurtler (see also the discussion in [0]). 



4.2 Angular Eigenstates and Monogenic Functions 

Returning to the Hamiltonian of equation H4.5|l , let us now consider the problem 
of a central potential V = V{r), ^ = 0, where r = \x\. We seek a set of angular- 
momentum operators which commute with this Hamiltonian. Starting with the 
scalar operator B-{x/\'V), where S is a spatial bivector, we find that 

[B-{xAV),n] = [B-{xAV),-jV] 
= ji7B-{xAV) 

^-jB-V. (4.27) 

But, since i?-V = [S, V]/2 and B commutes with the rest oiTC, we can rearrange 
the commutator into 

[B-{xAV) - ^B,H]=0, (4.28) 
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which gives us the required operator. Since B-{xf\W) — B /2 is an anti-Hermitian 
operator, we define a set of Hermitian operators as 

JB=]{B-{xhV)-\B). (4.29) 

The extra term of i^B is the term that is conventionally viewed as defining 'spin- 
f /2'. However, the geometric algebra derivation shows that the result rests solely 
on the commutation properties of the B-{xA'V) and V operators. Furthermore, 
the factor of one-half required in the Jb operators would be present in a space 
of any dimension. It follows that the factor of one-half in (|4.29(l cannot have 
anything to do with representations of the 3-D rotation group. 

From the STA point of view, Jb is an operator-valued function of the bivec- 
tor B. In conventional quantum theory, however, we would view the angular- 
momentum operator as a vector with components 

J^=U + is, (4.30) 

where Ei = {j /2)eijkjj^k- The standard notation takes what should be viewed 
as the sum of a scalar operator and a bivector, and forces it to look like the sum 
of two vector operators! As well as being conceptually clearer, the STA approach 
is easier to compute with. For example, it is a simple matter to establish the 
commutation relation 

[Jsi, Jsa] = -jJBtxB2, (4.31) 

which forms the starting point for the representation theory of the angular- 
momentum operators. 



The Spherical Monogenics 

The key ingredients in the solution of the Dirac equation for problems with 
radial symmetry are the spherical monogenics. These are Pauli spinors (even 
elements of the Pauli algebra H2.28|l 'l which satisfy the eigenvalue equation 

-xAVt/j = lip. (4.32) 

Such functions are called spherical monogenics because they are obtained from 
the 'monogenic equation' 

V* = (4.33) 

by separating into r'"ip{6,ip). Equation H4.33|l generalises the concept of an 
analytic function to higher dimensions |13l ; T . 

To analyse the properties of equation 14.32(1 we first note that 

[Jb,xAV] = 0, (4.34) 

which is proved in the same manner as equation (|4.28|) . It follows that ip can si- 
multaneously be an eigenstate of the a;AV operator and one of the Jb operators. 
To simplify the notation we now define 

Jkip = Jicki^ = {i(7k-{xhV) - \i(Jk)il)i(J-i. (4.35) 

We choose -0 to be an eigenstate of J3 , and provisionally write 

-ccAVV' = JsV' = mV'- (4.36) 
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Before proceeding, we must introduce some notation for a spherical-polar 
coordinate system. We define the {r, 9, cj)} coordinates via 

r = \/x^, cos9 — (Tj-x/r, tancf) ^ a2-x/(Ti-x. (4.37) 

The associated coordinate frame is 

e-r = sin0(cos0cri + sine/) (72) + coshers 

eg = rcos6'(cos0(Ti + sin0(T2) — rsinOas (4.38) 
— r 8m9{~ sin0 ai + coscj) 02). 

From these we define the orthonormal vectors {ur, erg, cr^} by 

ae - ee/r (4.39) 
0"^ = e0/(r sin6'). 

The {ar, erg, cr^} form a right-handed set, since 

CrO-eCT^ = i. (4.40) 

The vector ar satisfies 

XAV<7r = 2ar. (4.41) 

It follows that 

-xAViaripas) = -{I + 2)aripa3 (4.42) 

so, without loss of generality, we can choose / to be positive and recover the 
negative-Z states through multiplying by ar- In addition, since 

xAV{xAVi;) = l^il) (4.43) 

we find that 

19/ d%b\ 1 9V 

as Hi) + 

Hence, with respect to a constant basis for the STA, the components of i> are 
spherical harmonics and I must be an integer for any physical solution. 

The next step is to introduce ladder operators to move between different 
J3 eigenstates. The required analysis is standard, and has been relegated to 
Appendix^ The conclusions are that, for each value of I, the allowed values of 
the eigenvalues of J3 range from {I + 1/2) to —{I + 1/2). The total degeneracy 
is therefore 2{l + 1). The states can therefore be labeled by two integers I and 
m such that 

-ajAVV-r = ^V-r I > (4.45) 
JaV-i™ = (m + i)?A,™ -l-l<m<l. (4.46) 

Labelling the states in this manner is unconventional, but provides for many 
simplifications in describing the properties of the "0™. 

To find an explicit expression for the ipY^ we start from the highest-m eigen- 
state, which is given by 

V-l = sin' 61 (4.47) 
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I 


Eigenvalues of J3 


Degeneracy 


2 


5/2 5/2 


6 


1 


3/2 3/2 


4 





1/2 1/2 


2 


(-1) 


? 


7 


-2 


1/2 1/2 


2 



Table 5: Eigenvalues and degeneracies for the ■0™ monogenics. 



and act on this with the lowering operator J_. This procedure is described in 
detail in Appendix^ The result is the following, remarkably compact formula: 

- [{I + m+ l)P™(cos6l) - P;"+i(cos6i)i(T^]e"''^^'^3^ (4.48) 

where the associated Legendre polynomials follow the conventions of Grad- 
shteyn & Ryzhik 01]. The expression (|4.48|l offers a considerable improvement 
over formulae found elsewhere in terms of both compactness and ease of use. 
The formula H4.48|l is valid for non-negative I and both signs of to. The positive 
and negative m-states are related by 



^r(— .) = (4.49) 



The negative-Z states are constructed using H4.42II and the J3 eigenvalues are 
unchanged by this construction. The possible eigenvalues and degeneracies are 
summarised in TableEl One curious feature of this table is that we appear to be 
missing a line for the eigenvalue I = —1. In fact solutions for this case do exist, 
but they contain singularities which render them unnormalisable. For example, 
the functions 

have I = —1 and J3 eigenvalues +1/2 and —1/2 respectively. Both solutions are 
singular along the z-axis, however, so are of limited physical interest. 

4.3 Applications 

Having established the properties of the spherical monogenics, we can proceed 
quickly to the solution of various problems. We have chosen to consider two — 
the standard case of the Hydrogen atom, and the 'Dirac Oscillator' (45| . 

The Coulomb Problem 

The Hamiltonian for this problem is 

Zo. - 

n')p=P'ip ip + mip, (4.51) 

r 
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where a = e^/47r is the fine-structure constant and Z is the atomic charge. 
Since the Jb operators commute with 7i, -0 can be placed in an eigenstate of 
J3. The operator JiJi must also commute with H, but V does not, so both 
the V'/" and (TrtpY^(73 monogenics are needed in the solution. 
Though a;A V does not commute with H, the operator 

X = 7o(l-a;AV) (4.52) 

does, as follows from 

[70(1 - ccA V), V] = 270V - 7oViA V 

= 0. (4.53) 

We can therefore work with eigenstates of the K operator, which means that 
the spatial part of ^ goes either as 

ij{x, l + l)^ i^r^ir) + artpr^ir^icrs (4.54) 



or as 

ijix, -{I + 1)) = ari^r^Mr) + i^rHr)- (4.55) 

In both cases the second label in ip{x, I + 1) specifies the eigenvalue of K. The 
functions u{r) and v{r) are initially 'complex' superpositions of a scalar and an 
ias term. It turns out, however, that the scalar and ia^ equations decouple, 
and it is sufficient to treat u(r) and v{r) as scalars. 

We now insert the trial functions (|4.54|) and (|4.55|) into the Hamiltonian of 
equation H4.51|l . Using the results that 

-wr = i/rar^pr, -va.vr = -a + 2)/r vr, (4.56) 

and looking for stationary-state solutions of energy E, the radial equations 
reduce to 



u'\ _ f {K-l)/r -{E + Za/r + m)\ fu 
v' J \E + Za/r — m (— k— l)/r J \v 



(4.57) 



where k is the eigenvalue of K . {k is a non-zero positive or negative integer.) 
The solution of these radial equations can be found in many textbooks (see, 
for example, [351 137[l^ ). The solutions can be given in terms of confluent hy- 
pergeometric functions, and the energy spectrum is obtained from the equation 



E^ = m 



<^'"' (4.58) 



n2 + 2ni^ + (/ + 1)2 
where n is a positive integer and 

i^=[{l + lf + {Zaf]^/^. (4.59) 

Whilst this analysis does not offer any new results, it should demonstrate 
how easy it is to manipulate expressions involving the spherical monogenics. 
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The Dirac 'Oscillator' 

The equation describing a Dirac 'oscillator' was introduced as by Moshinsky 
and Szczepaniak as recently as 1989 |^|- The equation is one for a chargeless 
particle with an anomalous magnetic moment which, in the STA, takes the 
form 

Vi/'jo'a — ifiFipj^ — mtpjQ. (4.60) 

This equation will be met again in Section |51 where it is used to analyse the ef- 
fects of a Stern-Gerlach apparatus. The situation describing the Dirac oscillator 
is one where the F field exerts a linear, confining force described by 

F = X. (4.61) 

The Hamiltonian in this case is 

Hi/j = pV' ~ jmujx-ip + mil}. (4.62) 

It is a simple matter to verify that this Hamiltonian commutes with both the 
Jb and K operators defined above, so we can again take the wavefunction to 
be of the form of equations (|4.54|l and (j4.55|l . The resulting equations in this 

case are 



{K—l)/r — mujr --(E + m) \ f 

E^m [—K — l)/r + mujrJ \v 



(4.63) 



The equations arc simplified by transforming to the dimcnsionless variable p, 

p = {mcjy/'^r, (4.64) 

and removing the asymptotic behaviour via 

u = p'e^^'/^Mi (4.65) 

V ^ p^e^P^/'^U2. (4.66) 

The analysis is now slightly different for the positive- and negative-K equations, 
which we consider in turn. 

Positive-/^. The equations reduce to 



d f ui\ ( -{E + ■m)/^/mLu\ fui 

dp ^2/ ~ m)/V^ -2(; + l)/p + 2p) \U2 

which are solved with the power series 



(4.67) 



n=0 



"2 = }_J i?np'"+^ (4.69) 
n=0 



The recursion relations are 



2nAn^-^^^Br,-i (4.70) 
{2n + 2l + 3)Bn^^^^An + 2Bn-i, (4.71) 
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and the requirement that the series terminate produces the eigenvalue spectrum 



E-m=4nmuj n = l,2... (4.72) 

Remarkably, the energy levels do not depend on Z, so are infinitely degenerate! 
Negative-/^. In this case the equations reduce to 



dp \U2/ \{E - m) / y/muj 2p J \ U2 



(4.73) 



and are solved with the power series 

The recursion relations become 

{2n + 2l + 3)A^^^^^^Bn (4.76) 



2nBn = -^=An-i + 2S„_i, (4.77) 



and this time the eigenvalues are given by the formula 

E'^-m^ =2(2n + 2/ + l)rncj n = l,2.... (4.78) 

The energy spectrum only contains E through equations for E^. It follows 
that both positive and negative energies are allowed. The lowest positive-energy 
state (the ground state) has E^ — m'^ + Amuj, leading to a non-relativistic energy 
of ~ 2hu!. The groundstate is infinitely degenerate, whereas the first excited 
state has a degeneracy of two. The energy spectrum is clearly quite bizarre 
and does not correspond to any sensible physical system. In particular, this 
system does not reduce to a simple harmonic oscillator in the non-relativistic 
limit. The simple, if bizarre, nature of the energy spectrum is obscured in other 
approaches ^HlEZli which choose a less clear labeling system for the eigenstates. 



5 Propagators and Scattering Theory 

In this section we give a brief review of how problems requiring propagators 
are formulated and solved in the STA. The STA permits a first-order form of 
both Maxwell and Dirac theories involving the same differential operator — the 
vector derivative V. The key problem is to find Green's functions for the vector 
derivative that allow us to propagate initial data off some surface. We therefore 
start by studying the characteristic surfaces of the vector derivative. We then 
turn to the use of spinor potentials, which were dealt with in greater detail in 
The section concludes with a look at single-particle scattering theory. Using the 
mappings established in Sections 13 and 21 it is a simple matter to reformulate 
in the STA the standard matrix approach to scattering problems as described 
in |35l I37| . The STA approach allows for a number of improvements, however, 
particularly in the treatment of spin. This is a subject which was first addressed 
by Hestenes ^Hl; and our presentation closely follows his work. 
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5.1 Propagation and Characteristic Surfaces 



One of the simplest demonstrations of the insights provided by the STA formu- 
lation of both Maxwell and Dirac theories is in the treatment of characteristic 
surfaces. In the second-order theory, characteristic surfaces are usually found 
by algebraic methods. Here we show how the same results can be found using a 
simple geometric argument applied to first-order equations. Suppose, initially, 
that we have a generic equation of the type 

V^ = /(V',x), (5.1) 

where ip{x) is any multivector field (not necessarily a spinor field) and f{ip, x) 
is some arbitrary, known function. If we are given initial data over some 3-D 
surface are there any obstructions to us propagating this information off the 
surface? If so, the surface is a characteristic surface. We start at a point on 
the surface and pick three independent vectors {a, b, c} tangent to the surface 
at the chosen point. Knowledge of ip on the surface enables us to calculate 

a-V?/;, b-Vi/j and c-V^j. (5.2) 

We next form the trivector a/\b/\c and dualise to define 

n = iahbhc. (5.3) 

We can now multiply equation H5.1() by n and use 

= n-Vtp + i{aAbAc)-V^p (5.4) 
n-Wtp + {{aAbc-Wtp - aAcb-W^J + bAca-Wtp), (5.5) 

to obtain 



n-VV' = nf{'i{j,x) - i{aAbc-\/ip - ahcb-Vip + bAca-\/ip). (5.6) 

All of the terms on the right-hand side of equation l|5.6(l are known, so we can 
find n-W4> and use this to propagate ■0 in the n direction {i.e. off the surface). 
The only situation in which we fail to propagate, therefore, is when n remains 
in the surface. This occurs when 

nA{aAbAc) = 
=> nA{ni) = 

=^ n-n^O. (5.7) 

Hence we only fail to propagate when = 0, and it follows immediately that the 
characteristic surfaces of equation H5.1|l are null surfaces. This result applies to 
any first-order equation based on the vector derivative V, including the Maxwell 
and Dirac equations. The fundamental significance of null directions in these 
theories is transparent in their STA form. Furthermore, the technique extends 
immediately to a gravitational background, as described in j4Hj . 
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5.2 Spinor Potentials and Propagators 

A simple method to generate propagators for the Dirac theory is to introduce 
a spinor potential satisfying a scalar second-order equation. Suppose that ^ 
satisfies the Dirac equation 

Vipio-^ — rmpjo = 0. (5.8) 
ip can be generated from the (odd multivector) potential via 

Tp = Vcftia^, + m07o (5.9) 

provided that 

(V2+m2)0 = O. (5.10) 

The standard second-order theory can then be applied to (f>, and then used to 
recover "0. In |S] this technique was applied to constant-energy waves 

^ ^ i;{x)e~'''^^K (5.11) 
The Dirac equation then becomes 

V^pia3 + Ei;-m4!^0 (5.12) 

which is solved by 

ip = -V(l)ia3 + E(j) + m4> (5.13) 

where 

1 r pi(T3pr 

(bix) = -— (P\dS'\n'^P{x')- . (5.14) 

471 J r 

In this integral the initial data ip{x') is given over some closed spatial surface 
with normal n' = n{x'), and p and r are defined by 



p = \/ + HI? and r=|a; — a;'|. (5.15) 

Similar techniques can be applied to the propagation of electromagnetic waves 
(see for details). 



5.3 Scattering Theory 

We finish this short section with a brief look at how the matrix approach to scat- 
tering theory is handled in the STA, closely following the work of Hcstenes |48| . 
We continue to employ the symbol j for ia^ in places where it simplifies the 
notation. In particular, we employ the j symbol in the exponential terms intro- 
duced by Fourier transforming to momentum space. Where the j's play a more 
significant geometric role they are left in the ia^ form. 

We start by replacing the Dirac equation (|3.46|l with the integral equation 

ip{x) = V^(a;) + e j dV Sf[x - x')A{x')^{x') (5.16) 

where V'i is the asymptotic in-state which solves the free-particle equation, and 
Sf{x — x') is the STA form of the Feynman propagator. Substituting (|5.16fl 
into the Dirac equation, we find that Sf{x — x') must satisfy 

\I:rSF{x - x')M{x')ic73 - mSF(x - a;')M(a;')7o = ^{^ - x')M{x') (5.17) 
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for an arbitrary multivector M{x'). The solution to this equation is 

J (ztt)* — 

where, for causal propagation, the dE integral must arrange that positive- 
frequency waves propagate into the future {t > t') and negative-frequency waves 
propagate into the past [t' > t) . The result of performing the dE integral is 

SF{x-x')M = -e{t-t') I ^-l-(pM + mM7o)*a3e-J>-(^-^') 

J {It:)-^ IE 

+ 0{t' -t)j ^^{pM - mM7o)^<73eJ^'■(- " (5.19) 



where E = +y/p'^ + m? and M — M{x'). 
With ?/'diff(a;) defined by 

= ij{x) - ^,{x) (5.20) 
we find that, as t tends to +00, '0diff(2;) is given by 

V'diff'(x) = -e jd^x' j ^-L[pA(x')^(x')+mA(x')V'(x')7o]*a3e-J>-(^-^'). 

(5.21) 

We therefore define a set of final states ipf{x) by 

^f{x) = -e/ ^ [pfAix')^P{x') + mA{x')^P{x')^o]^cT3e~m-{^ - ^'), (5.22) 

which are plane- wave solutions to the free-field equations with momentum pf. 
ipdisix) can now be expressed as a superposition of these plane- wave states, 

d^Pf 
{2nf 



V'diff(^)= / T^^fi^)- (5.23) 



The Born Approximation and Coulomb Scattering 

In order to find ipf{x) we must evaluate the integral (|5.22() . In the Born ap- 
proximation, we simplify the problem by approximating ^{x') by Tpi{x'). In this 
case, since 

■<p^{x') ^iP,e~^P''''\ and TO?/'i7o = P^V'i, (5.24) 

we can write 

i'fix) ^~ej^ [pfA{x') + A(x')p#.»a3eJ'''-^'e-J>/-^, (5.25) 

where 

q^Pf-Pf (5.26) 

The integral in H5.25(l can now be evaluated for any given A-field. 

As a simple application consider Coulomb scattering, for which A{x') is given 

by 

Aix') = -^70. (5.27) 
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Inserting this in (|5.25|) and carrying out the integrals, we obtain 



^2 



i^fix) = ~Sf,i:^ia3^^S{Ef - E,) (5.28) 



E 

where 



Sfz = -;— [p/7o + loPt] / d a 



Z 
Q 



X ■ 



-^{2E + q). (5.29) 



Here E = Ef = Ei and a — /{An) is the fine-structure constant. The quantity 
Sfi contains all the information about the scattering process. Its magnitude 
determines the cross-section via 0S] 

'^^ Sf.Sfi (5.30) 



dflf 

and the remainder of Sfi determines the change of momentum and spin vectors. 
This is clear from (|5.28|l . which shows that Sfi must contain the rotor RfRi, 
where Ri and Rf are the rotors for the initial and final plane- wave states. 

Substituting (|5.29|l into H5.3()|l we immediately recover the Mott scattering 
cross-section 

— = ^{'iE^ - q^) = ^ (1 - sin2(^/2)) , (5.31) 

where 

= {Pf - Ptf = V(l - cos6l) and (3 = \p\/E. (5.32) 

The notable feature of this derivation is that no spin sums are required. Instead, 
all the spin dependence is contained in the directional information in Sfi. As well 
as being computationally more efhcient, the STA method for organising cross- 
section calculations offers deeper insights into the structure of the theory. For 
example, for Mott-scattering the directional information is contained entirely in 
the quantity 0H1 

Sf^ = - [p/70 + 7oP^] = L) + (5.33) 
■'to ^ 

where Lf and Li are the boosts contained in i?/ and Ri respectively. The 
algebraic structure 

Sf,=pfM + Mp^, (5.34) 

where M is some odd multivector, is common to many scattering problems. 

Since Sfi contains the quantity RfRi, we obtain a spatial rotor by removing 
the two boost terms. We therefore define the (unnormalised) rotor 

U'f, = LfiLj + Lf)L, = LfU + LfL„ (5.35) 

so that Ufi = U'j^/\Ufi\ determines the rotation from initial to final rest-spins. 
A simple calculation gives 

U'f,^2[{E + mf+PfPl (5.36) 
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Figure 1: Plane waves at a potential step. The spatial component of the 
wavevector lies in the x-y plane and the step lies in the y-z plane. 



hence the rest-spin vector precesses in the Pf/\Pi plane through an angle 5, 
where 

sin ^ 

tan(<5/2) = 77^ 7777^ ^ ^. (5.37) 

^' ' {E + m)/{E -m) +cosd ^ ' 

Whilst the derivations of the Mott scattering formula and the polarisation 
precession angle are only presented in outline here (further details are con- 
tained in ^1H|) it should be clear that they offer many advantages over the usual 
derivations [SSIISZI- AH the features of the scattering are contained in the single 
multivector Sji, algebraic form of which is very simple. Much work remains, 
however, if these techniques are to be extended to the whole of QED. 



6 Plane Waves at Potential Steps 

We now turn to a discussion of the matching of Dirac plane waves at a potential 
step. The case of perpendicular incidence is a standard problem and is treated in 
most texts [23 OH BHI ■ In order to demonstrate the power of the STA approach 
we treat the more general case of oblique incidence, adapting an approach used 
in electrical engineering to analyse the propagation of electromagnetic waves. 
A number of applications are given as illustrations, including the tunnelling of 
monochromatic waves and spin precession on total reflection at a barrier. We 
conclude the section with a discussion of the Klein paradox. 

The problem of interest is that of plane waves incident on a succession of 
potential steps. The steps are taken as lying along the x direction, with infinite 



38 



extent in the y and z directions. Since the spatial components of the incoming 
and outgoing wavevectors he in a single plane, the matching problem can be 
reduced to one in two dimensions. The analysis is simplified further if the 
wavevectors are taken to lie in the ia^ plane. (Other configurations can always 
be obtained by applying a rotation.) The arrangement is illustrated in Figure^ 
The waves all oscillate at a single frequency E, and the Dirac equation in the 
ith region is 

(E - eVi)^! ^ -V^jiaz + mjoTpjo- (6.1) 

By continuity of ip at each boundary, the y component of the wavevector, py, 
must be the same in all regions. For the ith region we define 

E[ = E~ eV^ (6.2) 

and, depending on the magnitude of Vi, the waves in the this region will be 
either travelling or evanescent. For travelling waves we define (dropping the 
subscripts) 



pI = E''^ ~ pI - \E-eV\> ^Jpl + m?. (6.3) 

In terms of the angle of incidence (j) we also have 

Px = pcoscf), Py — psini/), E'"^ = p^ + rn? . (6.4) 
For evanescent waves we write 



K 



E"-+pl+m^ \E- eV\< Jpl+'m?. (6.5) 



In all the cases that we study, the incoming waves are assumed to be positive- 
energy travelling waves in a region where eV < E — \J + 'ti^ . Recalling the 
plane- wave solutions found in Section 13.31 the travelling waves are given by 

V-* = [cosh(w/2) + sinh(?i/2)(cos0(Ti +sin<?!>CT2)]$e~*<^3(^*~^-^^P«3')r 

(6.6) 

i)'' = [cosh(u/2)-f sinh(u/2)(-cos(/iCTi -Ksin(/)CT2)]$e-^^3(^* + P-^-^'«2')i?, 

(6.7) 

where 

tanh(w/2) = p/(E' + m) (6.8) 

so that 

sinhu = p/m, coshw = E' /m. (6.9) 

The transmission and reflection coefficients T and R are scalar -\-iaz combina- 
tions, always appearing on the right-hand side of the spinor. The fact that py 
is the same in all regions gives the electron equivalent of Snell's law, 

sinhu sin0 ~ constant. (6.10) 

The Pauli spinor $ describes the rest-spin of the particle, with <I> = 1 giving 
spin-up and $ = — icr2 spin down. Other situations are, of course, built from 
superpositions of these basis states. For these two spin basis states the spin 
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vector is ±173, which Ucs in the plane of the barrier and is perpendicular to 
the plane of motion. Choosing the states so that the spin is aligned in this 
manner simplifies the analysis, as the two spin states completely decouple. Many 
treatments (including one published by some of the present authors [H]) miss 
this simplification. 

There are three matching situations to consider, depending on whether the 
transmitted waves are travelling, evanescent or in the Klein region {eV > E + 

'\Jpy + rn^)- We consider each of these in turn. 

6.1 Matching Conditions for Travelling Waves 

The situation of interest here is when there are waves of type (lti.7|l in both 
regions. The matching condition in all of these problems is simply that ip 
is continuous at the boundary. The work involved is therefore, in principle, 
less than for the equivalent non-relativistic problem. The matching is slightly 
different for the two spins, so we consider each in turn. 

Spin-Up ($ = 1) 

We simplify the problem initially by taking the boundary at a; = 0. Steps at 
other values of x are then dealt with by inserting suitable phase factors. The 
matching condition at x = reduces to 

[cosh(ui/2) + smh{u j2)aie'l''''^^Tj + [cosh(M,/2) - sinh(u,/2)CTie~'^»*'^«]i?J 

= [cosh(7.,+i/2) + sinh(«,+i/2)aie'^«+i^'^3]2;T^^ 

+ [cosh(u,+i/2) - sinh(u,+i/2)(7ie-'^«+i*'^3]i?J+i. 

Since the equations for the refiection and transmission coefficients involve only 
scalar and ias terms, it is again convenient to replace the ia^ bivector with the 
symbol j. If we now define the 2x2 matrix 

/ cosh(Mi/2) cosh(ui/2) 
' " Unh(M»/2)eJ>' - sinh(w,/2)e~^>' 

we find that equation H6.11|l can be written concisely as 

The Ai matrix has a straightforward inverse, so equation 16.13|l can be easily 
manipulated to describe various physical situations. For example, consider plane 
waves incident on a single step. The equation describing this configuration is 



(6.12) 



simply 



^1 ri) = ^2 (^f] (6.14) 



so that 



R[J sinhwi cos(/)i 

/ sinh(ui/2) cosh(M2/2)e-J''^i + cosh(Mi/2) sinh(M2/2)e-?''^2' 
I sinh(ui/2)cosh(M2/2)e-?''^i - cosh(ui/2) sinh(M2/2)e-?''^2 



(6.15) 
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from which the reflection and transmission coefficients can be read off. The 
case of perpendicular incidence is particularly simple as equation H6.13|l can be 
replaced by 



Ti+i\ _ 1 / sinhi(ui+i + u^) sinhi(ui+i - Ui)\ f Ti 
Ri+iJ sinhuj+i Vsinh|(iii+i - Wj) sinh2 (li^+i + l^i?, 



(6.16) 



which is valid for all spin orientations. So, for perpendicular incidence, the 
reflection coefficient r — and transmission coefficient t — T2/T1 at a 

single step are 

sinhi(ui -Us) sinhui 
r = f- -, t = T- -, (d-1') 

sinh2(ui + U2) sinh2(ui + U2) 

which agree with the results given in standard texts (and also |36p. 
Spin-Down ($ = — io'2) 

The matching equations for the case of opposite spin are 

[cosh(ui/2) + sinh(ui/2)crie''^"*'^3](-icr2)T;^ 
+ [cosh(ui/2) - smh{uj2)aie-'^'^'^'^']{-ia2)Ri 

= [cosh(u,+i/2) + sinh(u,+i/2)aie<^>+i»'^3](_i^2)T;Vi 

+ [cosh(u,+i/2) - sinh(u,+i/2)aie-'^'+i^'^3](-ifT2)ii'j+i. (6.18) 

Pulling the ia2 out on the right-hand side just has the effect of complex- 
conjugating the reflection and transmission coefficients, so the matrix equa- 
tion H6.11|l is unchanged except that it now relates the complex conjugates of 
the reflection and transmission coefficients. The analog of equation (|6.13(l is 
therefore 

As mentioned earlier, the choice of alignment of spin basis states ensures that 
there is no coupling between them. 

One can string together series of barriers by including suitable 'propagation' 
matrices. For example, consider the set-up described in Figure|21 The matching 
equations for spin-up are, at the first barrier. 



J,, J -A, (6.20) 

and at the second barrier. 
Equation H6.21|l demonstrates neatly how matrices of the type 

(6.22) 
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eV 




^ 1 X 

X = X = d 

Figure 2: Plane waves scattering from a barrier. The barrier has height eV and 
width d. Quantities inside the barrier are labeled with a subscript 1, and the free 
quantities have no subscripts. The phases are given by 5i = md smhui cos(j)i 
and S = Pxd. 



where S — Pxd and d is the distance between steps, can be used to propagate 
from one step to the next. In this case the problem is reduced to the equation 



1 

-1 







coaSiI — j s\n5i A ^Ai 
which quickly yields the reflection and transmission coefflcients. 



(6.23) 



6.2 Matching onto Evanescent Waves 

Before studying matching onto evanescent waves, we must first solve the Dirac 
equation in the evanescent region. Again, the two spin orientations behave 
differently and are treated separately. Taking spin-up first and looking at the 
transmitted (decaying) wave in the evanescent region, the solution takes the 
form 

= [cosh(M/2) + sinh(M/2)a2]e-''^e-^'^3(£t - Pyy)^^ (6.24) 
Substituting this into the Dirac equation yields 

(i?'7o + Pvl2)e^'''l'^ = e^'^^l'^im^^ - n^^) (6.25) 

which is consistent with the definition of n 16.5|l . From equation H6.25|l we find 
that 

which completes the solution. For the incoming (growing) wave we flip the sign 
of K. We therefore define via 

tanh(«V2) - 1^ (6.27) 

hi.' _L rn 
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and write the outgoing and incoming spin-up waves in the evanescent region as 
V-* = [cosh(u+/2) + sinh(u+/2)CT2]e"''^e-*'^3(i?t-P,j/)7^T (6.28) 
V^'' = [cosh(w-/2) + sinh(u-/2)a2]e''^e"*'^3(^* " Pyy^ R\ (6.29) 
If we now consider matching at x = 0, the continuity equation becomes 

[cosh(u,/2) + sinh(u,/2)crie'?^»^'^3]7:T 
+ [cosh(u,/2) - sinh(w,/2)crie-'?^»*'^3]i^T 

= [cosh(u+ i/2) + sinh(ii+ i/2)a2]i;Vi 
+ [coshKr^i/2) + sinhKr^i/2)a2)]i?J+i. 

On defining the matrix 



+ _ / cosh(w+/2) cosh(w^ /2) , 

~ ^^jsinh(u+/2) +jsinh(ur/2)^) ' ^'^•'^^^ 



BU ) ■ (6-32) 



we can write equation (|6.3UI) compactly as 



Again, either of the matrices can be inverted to analyse various physical situa- 
tions. For example, the case of total reflection by a step is handled by 

from which one finds the reflection coefficient 

_ tanh(M+/2) + tanh(M/2)jeJ''^ 
^ tanh(u+/2) -tanh(w/2)je-J> 

which has | = I, as expected. The subscripts on ui, and 0i are all obvious, 
and have been dropped. 

The case of spin-down requires some sign changes. The spinors in the evanes- 
cent region are now given by 

ij/ = [cosh(u-/2) +sinh(u-/2)c72](-«CT2)e~''^e~^^3(i5i-P,y)yi (6.35) 
= [cosh(w+/2) + sinh(ii+/2)(72](-ia2)e''^e-^'^3(£t-p«y)^i (6.36) 

and, on defining 

_ / cosh(u,-/2) cosh(w+/2) \ , , 

^ " l^jsinh(ur/2) -j sinh(M+/2)y' ' ^^'^'^ 

the analog of equation H6.32|l is 

These formulae are now applied to two situations of physical interest. 
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6.3 Spin Precession at a Barrier 



When a monochromatic wave is incident on a single step of sufficient height that 
the wave cannot propagate there is total reflection. In the preceding section we 
found that the reflection coefficient for spin-up is given by equation (|6.34(l , and 
the analogous calculation for spin-down yields 

^1 _ tanh(u-/2) - tanh(u/2)je~J> 
^ ta.-Dh{u- /2) + ta.nh{u/2)je3<t> 

Both and are pure phases, but there is an overall phase difference between 
the two. If the rest-spin vector s = <I>cr3<I> is not perpendicular to the plane of 
incidence, then this phase difference produces a precession of the spin vector. 
To see how, suppose that the incident wave contains an arbitrary superposition 
of spin-up and spin-down states, 

$ = cos(0/2)e^'^3<^i - sin(^/2)ia2e*'^3'^2 = ^i<J-i4>/2^-ia2e/2^-i<j-ie/2 ^ (g 
where 

cf)^ cf)i- cf)2, e = (t)i+(f>2, (6.41) 

and the final pure-phase term is irrelevant. After reflection, suppose that the 
separate up and down states receive phase shifts of (5^ and respectively. The 
Pauli spinor in the reflected wave is therefore 

= cos(0/2)e^'^3(<^i + 5^) _ sin(0/2)z(T2e*'^3(02 + S^) 

^ ^i(J35 /2 ^iaz(t)/2 ^~ia29 /2 ^-iaze /2 ^ (6.42) 



where 



5 = 5^-5^ (6.43) 



and again there is an irrelevant overall phase. The rest-spin vector for the 
reflected wave is therefore 

so the spin- vector precesses in the plane of incidence through an angle (5^ — 5-^ . 
If (5^ and 5-^ are deflned for the asymptotic (free) states then this result for the 
spin precession is general. 

To find the precession angle for the case of a single step we return to the 
formulae (|6.34|l and H6.39|l and write 

^3 ^ — ^ T 

(tanh(M+/2) + tanh(u/2)je^'<?^)(tanh(u-/2) + tanh(u/2)je^''^) 

~ (tanh(u+/2) - tanh(u/2)je":''':*)(tanh(u-/2) - ta.nh{u/2)je~i't>) 

(6.45) 

If we now recall that 

tanh(V2) = pI{E + m), tanh(^) = ^^"^^^^ , (6.46) 

2 hj — eV + m 
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we find that 

m COS0 + jE svncj) 

The remarkable feature of this result is that all dependence on the height of the 
barrier has vanished, so that the precession angle is determined solely by the 
incident energy and direction. To proceed we write 

TO cos4> — jE sin0 = pe^'^ (6.48) 

so that 

tana = — coshu tani/). (6.49) 

Equation H6.47|l now yields 

tan((5/2 — 0) = — coshw taiKp, (6.50) 

from which we obtain the final result that 

= _ (cosh.-l)tan0 _ 
^ ' ^ l + cosh-utan2(/) ^ ' 

A similar result for the precession angle of the rest-spin vector was obtained 
by Fradkin & Kashuba [101 using standard techniques. Readers are invited to 
compare their derivation with the present approach. The formula H6.51|l agrees 
with equation (|5.37(l from Section [E31 since the angle 6 employed there is related 
to the angle of incidence </> by 

6l==7r-20. (6.52) 

Since the decomposition of the plane-wave spinor into a boost term and a 
Pauli spinor term is unique to the STA, it is not at all clear how the conventional 
approach can formulate the idea of the rest-spin. In fact, the rest-spin vector is 
contained in the standard approach in the form of the 'polarisation operator' j5()l 
I3T] which, in the STA, is given by 

j 

0{n) = y[ipp-n + ■yoi{n/\p)-p] (6.53) 

P 

where n is a unit spatial vector. This operator is Hermitian, squares to 1 and 
commutes with the free-field Hamiltonian. If we consider a free-particle plane- 
wave state, then the expectation value of the 0{ui) operator is 

™ {^L{a.■ppL<^a, + a.^ppL<,a,)) 



(iP^) Ep^ 



TO 

Ej^ 



[L{a^-ppL + a^AppL)]-s (6.54) 



where L is the boost 



E + m + p 
^2m{E + to) 

and p is replaced by its eigenvalue p. To manipulate equation (|6.54|) we use the 
facts that L commutes with p and satisfies 

l? = {E^p)lm (6.56) 
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X = —d X = Q 

Figure 3: Schematic representation of plane-wave tunnelling. 



-^X 



to construct 

TO 



[L{ai-ppL + atAppL)] 



m , (E + p) 

— Wi-P P CFiPP 



Ep 

E — m 
2m 



TO 



{E + m + p)a,{E + m - p)]. (6.57) 



Since only the relative vector part of this quantity is needed in equation (|6.54|l 
we are left with 



E-m 
2Ep^ 



{2ai pp + (E + mfdi - paip) = 



1 



2E{E + m) 



{p^a, + [E + mfai) 



(6.58) 



The expectation value of the 'polarisation' operators is therefore simply 



(6.59) 



which just picks out the components of the rest-spin vector, as claimed. 

For the case of the potential step, 0{n) still commutes with the full Hamilto- 
nian when n is perpendicular to the plane of incidence. In their paper Fradkin & 
Kashuba decompose the incident and reflected waves into eigenstates of 0{n), 
which is equivalent to aligning the spin in the manner adopted in this section. 
As we have stressed, removing the boost and working directly with <I> simpli- 
fies many of these manipulations, and removes any need for the polarisation 
operator. 



6.4 Tunnelling of Plane Waves 

Suppose now that a continuous beam of plane waves is incident on a potential 
barrier of finite width. We know that, quantum-mechanically, some fraction of 
the wave tunnels through to the other side. In Section[71we address the question 
'how long does the tunnelling process take?'. To answer this we will need to 
combine plane- wave solutions to construct a wavepacket, so here we give the 
results for plane waves. The physical set-up is illustrated in Figure El 
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The matching equation at the x — boundary is, for spin-up, 







(6.60) 



where the A and B^ matrices are as defined in H6.12|l and (|6.31|l respectively. 
All subscripts can be dropped again as A always refers to free space and B^ to 
the barrier region. The matching conditions at a; = —d require the inclusion of 
suitable propagators and the resulting equation is 









(6.61) 



Equation (|6.61|l shows that the relevant propagator matrix for evanescent waves 
is 



e-d A 
e-^'^j ■ 

The problem now reduces to the matrix equation 



(6.62) 



cosh(Kd) + sinh(Kd)A ^B^ 



-1 



B^ 







(6.63) 

from which the reflection and transmission coefficients are easily obtained. 

The applications to tunnelling discussed in Section |7| deal mainly with per- 
pendicular incidence, so we now specialise to this situation. For perpendicular 



incidence we can set u 



+ 



u' where 



tanhii' 



E' 



(6.64) 



It follows that the equations for spin-up and spin-down are the same, and we 
can remove the up-arrows from the preceding equations. Equation (|6.60(l now 
yields 



Ri 



T /sinh(u72) cosh(M/2) - j cosh(M72) sinh(w/2) 
sinhu' I sinh(u72) cosh(u/2) -I- j cosh(u72) sinh(w/2) 



(6.65) 



and from Ti and i?i the current in the evanescent region can be constructed. 
The ratio Ji/Jo may be interpreted as defining a 'velocity' inside the barrier. 
The consequences of this idea were discussed in jSj where it was concluded that 
the tunnelling times predicted by this velocity are not related to tunnelling 
times measured for individual particles. The reasons for this are discussed in 
Section 13 

Multiplying out the matrices in equation H6.63|l is straightforward, and yields 




= T' 



cosh(K(i) — j sinh(K(i)(£'£" — w?)/{kp) 



which solves the problem. The transmission coefficient is 

-jdp 



t 



Kpe 



Kpcosh(K(i) — jijp — eVE) smhi^nd) 
which recovers the familiar non-relativistic formula in the limit E 



(6.66) 



(6.67) 



m. 
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6.5 The Klein Paradox 



In the Klein region, eV — E > yPy + travelling wave solutions exist again. 
To find these we observe that plane-wave solutions must now satisfy 

(p — 6^70)^' — mipjo (6.68) 

and, as p — eV"/Q has a negative time component, ipip must now be —1. We could 
achieve this flip by inserting a '/3-factor' of the type described in Section H3.3|l . 
but this would mix the rest-spin states. It is more convenient to work with 
solutions given by 

V'* = [cosh(V2) + sinh(u/2)(cos0(Ti - sin0(72)]ai$e-^'^3(i5t + P.a; - 

(6.69) 

^P'' = [cosh(V2) + sinh(u/2)(-cos(/.ai - sin0cr2)]c7i$e-*^3(^* " " ^^S^^i?, 

(6.70) 

where the choice of cti or a2 on the right-hand side of the boost is merely a 
phase choice. To verify that tp* is a solution we write the Dirac equation as 

[{E - eVho - P.li + P,72]e(^°^'^'^^ " ^^^^'^^)^/^a^ 
which holds provided that 

tanh(u/2) = \ (6.72) 
m + eV — E 

It follows that 

m coshu = eT^ — i?, msinhw=p. (6.73) 

The current obtained from -0* is found to be 

V'*7oV'* = [aV - E)-/o + Pxli - Pvl2 (6.74) 

which is future pointing (as it must be) and points in the positive-x direction. 
It is in order to obtain the correct direction for the current that the sign of Px is 
changed in H6.69|l and H6.70|l . As has been pointed out by various authors 
1521 153| . some texts on quantum theory miss this argument and match onto 
a solution inside the barrier with an incoming group velocity |35l I37j . The 
result is a reflection coefficient greater than 1. This is interpreted as evidence 
for pair production, though in fact the effect is due to the choice of boundary 
conditions. To find the correct reflection and transmission coefficients for an 
outgoing current, we return to the matching equation which, for spin-up, gives 

[cosh(ui/2) + sinh(ui/2)CTie'^'^'^3]jTT 
+ [cosh(ui/2) - sinh(M,/2)crie-'?^"*'^3]^T 



= [cosh(u,+i/2) + sinh(u,+i/2)aie-<^'+i»'^3]^^2:T^^ 
+ [cosh(u,+i/2) - sinh(u,+i/2)aie'^'+i*'^3]^^i^T_^^. 

This time we define the matrix 



(6.75) 



fsmhiu,/2)ej'l'' - sinh(V2)e~^''^' , 
' \ cosh(M,/2) cosh(u,/2) ' ^ ' 
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so that equation l|t).75|l becomes 



..J\]=C,+Jlr]. (6.77) 



It should be noted that 

I 

The corresponding equation for spin-down is simply 



C.-f? J) A,. (6.78) 



i I = ) • (6-79) 



The Klein 'paradox' occurs at a single step, for which the matching equation 



Inverting the Ai matrix yields 

'tI\ _ tJ / cosh(M + u')/2 

' sinhu COS0 Vsinh(u/2) sinh(w72)e2i</' - cosh(u/2) cosh(M72) 

(6.81) 

from which the reflection and transmission coefficients can be read off. (The 
primed quantities relate to the barrier region, as usual.) In particular, for 
perpendicular incidence, we recover 

^ ^ cosh(M - u')/2 ^ ^ sinhu 

'"^ cosh(u + u')/2' ^ cosh(M + u')/2 ^' ' 

as found in j36l . The reflection coefflcient is always < 1, as it must be from 
current conservation with these boundary conditions. But, although a reflection 
coefficient < 1 appears to ease the paradox, some difflculties remain. In partic- 
ular, the momentum vector inside the barrier points in an opposite direction to 
the current. 

A more complete understanding of the Klein barrier requires quantum field 
theory since, as the barrier height is > 2m, we expect pair creation to occur. 
An indication that this must be the case comes from an analysis of boson modes 
based on the Klein-Gordon equation. There one finds that superradiance (r > 1) 
does occur, which has to be interpreted in terms of particle production. For the 
fermion case the resulting picture is that electron-positron pairs are created and 
split apart, with the electrons travelling back out to the left and the positrons 
moving into the barrier region. If a single electron is incident on such a step 
then it is reflected and, according to the Pauli principle, the corresponding 
pair-production mode is suppressed. 

A complete analysis of the Klein barrier has been given by Manogue |52| 
to which readers are referred for further details. Manogue concludes that the 
fermion pair-production rate is given by 

r . /-^ 5^ y |Tf (6.83) 
J (27r)2 J 2tt sinhu ' ^ ^ 
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where the integrals run over the available modes in the Klein region, and the 
sum runs over the two spin states. This formula gives a production rate per unit 
time, per unit area, and applies to any shape of barrier. The integrals in (|6.83() 
are not easy to evaluate, but a useful expression can be obtained by assuming 
that the barrier height is only slightly greater than 2m, 

eV ^2m{l + e). (6.84) 

Then, for the case of a single step, we obtain a pair-production rate of 

7rm e 

r = ^, (6.85) 

to leading order in e. The dimensional term is which, for electrons, cor- 
responds to a rate of 10'*^ particles per second, per square meter. Such an 
enormous rate would clearly be difhcult to sustain in any physically-realistic 
situation! 

The results obtained in this section are summarised in Table 



7 Tunnelling Times 

In this Section we study tunnelling phenomena. We do so by setting up a 
wavepacket and examining its evolution as it impinges on a potential barrier. 
The packet splits into reflected and transmitted parts, and the streamlines of the 
conserved current show which parts of the initial packet end up being transmit- 
ted. The analysis can be used to obtain a distribution of arrival times at some 
fixed point on the far side of the barrier, which can be compared directly with 
experiment. The bulk of this section is concerned with packets in one spatial 
dimension, and compares our approach to other studies of the tunnelling-time 
problem. The section ends with a discussion of the complications introduced in 
attempting 2- or 3-dimensional simulations. 

The study of tunnelling neatly combines the solutions found in SectionElwith 
the views on operators and the interpretation of quantum mechanics expressed 
in Section |3 Tunnelling also provides a good illustration of how simple it is 
study electron physics via the Dirac theory once the STA is available. 

7.1 Wavepacket Tunnelling 

In Section 16.41 we studied tunnelling of a continuous plane wave through a po- 
tential barrier. It was found that the growing and decaying waves in the barrier 
region are given by equations H6.29(l for spin-up and H6.36(l for spin-down. Re- 
stricting to the case of perpendicular incidence, the amplitudes of the reflected 
and transmitted waves are given by equation (|6.65|1 . It follows that, for arbitrary 
spin, the wavefunction in the barrier region is 

= [cosh(u72)$ + sinh(u72)cr2cr3$cr3]e"''^e~*'^3'^*a-f 
[cosh(u72)$ - sinh(u72)(T2(T3$a3]e''^e-^'^3i5tc,* (7.1) 

where 

T' 

a = . , , [sinh(u72) cosh(u/2) - ias cosh(u72) sinh(u/2)] (7.2) 
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Travelling Waves 



v^* = 


[cosh(u/2) + sinli(u/2)aie^*'^3]$e-*^3(^* " " 


-PvV) 


tp'' = 


cosli(u/2) - sinh(u/2)aie-'^*^3]$e-*^3(^* + P-^ 


-PvV) 




tanh(u/2) =p/{E-eV + m) 





Evanescent Waves 



v^* = 


[cosh(u± /2) + sinh(u± /2)cr2]e"''^e-*<^3(£^t 


-Pyy)T 


tp'' = 


[cosh(ii=F/2) + sinh(ii=F/2)a2]e'^^e-^^3(^* " 


-Pyy)R 




(upper/lower signs = spin up/down) 






tanh(u±/2) = {py ± k)/{E - eV + m) 





Klein Waves 





cosh(u/2) + sinh(u/2)crie-'?^*'^3]cr^$e 


-ia3{Et + pxX 


-PvV) 


v^*^ = 


[cos1i(m/2) - sinh(u/2)crie'^*'^;^]cri$e- 


-i(T3{Et -pxX - 


-PvV) 




tanh(M/2) = p/{m + eV — 


E) 





Matching Matrices 



^ _ / cosh(u/2) cosh(M/2) ^ 

~ l^sinh(u/2)e-?''?^ - sinh(M/2)e--?>^ 
B+ - ( cosh(w+/2) cosh(w-/2) \ 
I j sinh(u+/2) +j sinh(?i~/2) j 
sinh(u/2)e-?> - sinh(w/2)e--?>^ 
\ cosh(u/2) cosh('u/2) , 

A* , B~ , C* for spin down. 



Propagators 
\^ e-^'^P- j ' e-'^'^ y 



Table 6: Summary of results for plane waves incident on a potential step. The 

waves travel in the x — y plane and the steps lie in the y — z plane. The matching 
matrices relate T and R on either side of a step. 
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and 



tanhK/2)^-^^ I , n^ = m'-E'\ (7.3) 

The current in the barrier region is 
IT'P 

ijjijo^Pi = J — ^-\m'^eVcosh(2Kx) + E' {p^ - Eev) 

+ pK^ai — mneV smh{2Kx){i(7i)-s]jo ('^•4) 
from which we can define a 'velocity' 

dx J pK^ 



dt J-7o m'^eV cosh{2 Kx) + E'{p'^ — Eev)' 



(7.5) 



In fact, the velocity (|7.5() does not lead to a sensible definition of a tunnelling 
time for an individual particle As we shall see shortly, an additional phe- 
nomenon underlies wavepacket tunnelling, leading to much shorter times than 
those predicted from (|7.5|l . To study wavepacket tunnelling it is useful, initially, 
to simplify to a one-dimensional problem. To achieve this we must eliminate the 
transverse current in (|7.4() by setting s = ±ai. This is equivalent to aligning 
the spin vector to point in the direction of motion. (In this case there is no 
distinction between the laboratory and comoving spin.) With <i> chosen so that 
s = (Ti it is a now a simple matter to superpose solutions at t = to construct 
a wavepacket centred to the left of the barrier and moving towards the barrier. 
The wavepacket at later times is then reassembled from the plane-wave states, 
whose time evolution is known. The density J'^ =Jo-J can then be plotted as 
a function of time and the result of such a simulation is illustrated in Figure 21 
The Dirac current J = iplo^' is conserved even in the presence of an elec- 
tromagnetic field. It follows that J defines a set of streamlines which never end 
or cross. Furthermore, the time-component of the current is positive-definite 
so the tangents to the streamlines are always future-pointing timelike vectors. 
According to the standard interpretation of quantum mechanics, J^{x,t) gives 
the probability density of locating a particle at position x at time t. But, 
considering a flux tube defined by adjacent streamlines, we find that 

p{to,xo)dxo = p{ti,xi)dxi (7.6) 

where {t^^xo) and (ti,xi) are connected by a streamline. It follows that the 
density J'^ flows along the streamlines without 'leaking' between them. So, in 
order to study the tunnelling process, we should follow the streamlines from 
the initial wavepacket through spacetime. A sample set of these streamlines 
is shown in Figure A signiflcant feature of this plot is that a continuously- 
distributed set of initial input conditions has given rise to a disjoint set of 
outcomes (whether or not a streamline passes through the barrier). Hence the 
deterministic evolution of the wavepacket alone is able to explain the discrete 
results expected in a quantum measurement, and all notions of wavcfunction 
collapse are avoided. This is of fundamental significance to the interpretation of 
quantum mechanics. Some consequences of this view for other areas of quantum 
measurement have been explored by Dewdney et al. |28j and Vigier et al. |27j . 
though their work was founded in the Bohmian interpretation of non-relativistic 



52 




Figure 4: Evolution of the density J° as a function of time. The initial packet 
is a Gaussian of width Ak = 0.04A~^ and energy 5eV. The barrier starts at the 
origin and has width 5 A and height lOeV. The top line shows the density profile 
at times —0.5 x lO^^^s and —0.1 x 10~^*s, and the bottom line shows times 
0.1 X 10" ^^s and 0.5 x 10" ^''s. In all plots the vertical scale to the right of the 
barrier is multiplied by lO"' to enhance the features of the small, transmitted 
packet. 

quantum mechanics. The results presented here are, of course, independent of 
any interpretation — we do not need the apparatus of Bohm/de Broghe theory in 
order to accept the vahdity of predictions obtained from the current streamhnes. 

The second key feature of the streamhne plot in FigureElis that it is only the 
streamlines starting near the front of the initial wavepacket that pass through 
the barrier. Relative to the centre of the packet, they therefore have a 'head 
start' in their arrival time at some chosen point on the far side of the barrier. 
Over the front part of the barrier, however, the streamlines slow down consid- 
erably, as can be seen by the change in their slope. These two effects, of picking 
out the front end of the packet and then slowing it down, compete against each 
other and it is not immediately obvious which dominates. To establish this, 
we return to Figure 01 and look at the positions of the wavepacket peaks. At 
t = 0.5 X lO^^'^s, the peak of the transmitted packet lies at a; = 70A, whereas the 
peak of the initial packet would have been at x = 66A had the barrier not been 
present. In this case, therefore, the peak of the transmitted packet is slightly 
advanced, a phenomenon often interpreted as showing that tunnelling particles 
speed up, sometimes to velocities greater than c |54) . The plots presented here 
show that such an interpretation is completely mistaken. There is no speeding 
up, as all that happens is that it is only the streamlines from the front of the 
wavepacket that cross the barrier (slowing down in the process) and these re- 
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STREAMLINES 




time {10«t-14 sees) 

Figure 5: Particle streamlines for the packet evolution shown in Figure^ Only 
the streamlines from the very front of the packet cross the barrier, with the 
individual streamlines slowing down as they pass through. 

assemble to form a localised packet on the far side. The reason that tunnelling 
particles may be transmitted faster than free particles is due entirely to the 
spread of the initial wavepacket. 

There is considerable interest in the theoretical description of tunnelling 
processes because it is now possible to obtain measurements of the times in- 
volved. The clearest experiments conducted to date have concerned photon 
tunnelling j55| . where an ingenious 2-photon interference technique is used to 
compare photons that pass through a barrier with photons that follow an unob- 
structed path. The discussion of the results of photon-tunnelling experiments 
usually emphasise packet reshaping, but miss the arguments about the stream- 
lines. Thus many articles concentrate on a comparison of the peaks of the in- 
cident and transmitted wavepackets and discuss whether the experiments show 
particles travelling at speeds > c . As we have seen, a full relativis- 

tic study of the streamlines followed by the electron probability density show 
clearly that no superluminal velocities are present. The same result is true for 
photons, as we will discuss elsewhere. 

Ever since the possibility of tunnelling was revealed by quantum theory, peo- 
ple have attempted to define how long the process takes. Reviews of the various 
different approaches to this problem have been given by Hauge & St0vneng |F7) 
and, more recently, by Landauer & Martin [^. Most approaches attempt to 
define a single tunnelling time for the process, rather than a distribution of 
possible outcomes as is the case here. Quite why one should believe that it is 
possible to define a single time in a probabilistic process such as tunnelling is 
unclear, but the view is still regularly expressed in the modern literature. A fur- 
ther flaw in many other approaches is that they attempt to define how long the 
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particle spent in the barrier region, with answers ranging from the implausible 
(zero time) to the utterly bizarre (imaginary time). From the streamline plot 
presented here, it is clearly possible to obtain a distribution of the times spent 
in the barrier for the tunnelling particles, and the answers will be relatively 
long as the particles slow down in the barrier. But such a distribution neglects 
the fact that the front of the packet is preferentially selected, and anyway does 
not appear to be accessible to direct experimental measurement. As the recent 
experiments show |55) . it is the arrival time at a point on the far side of the 
barrier that is measurable, and not the time spent in the barrier. 

7.2 2-Dimensional Simulations 

In the preceding section we simplified the problem in two ways: by assum- 
ing perpendicular incidence, and by aligning the spin-vector in the direction of 
motion. For other configurations more complicated two-dimensional or three- 
dimensional simulations are required. As well as the obvious numerical com- 
plications introduced there are some further difficulties. For the 1-D plots just 
shown, there was no difficulty in deciding which part of the wavepacket was 
transmitted and which was reflected, since the bifurcation point occurred at 
some fixed value of x. For 2-D or 3-D simulations, however, the split of the 
initial packet into transmitted and reflected parts occurs along a line or over 
a 2-D surface. Furthermore, this split is spin-dependent — if one constructs 
a moving wavepacket, one finds that the streamlines circulate around the spin 
axis |3()j . (A similar circulation phenomena is found in the ground state of the 
hydrogen atom [S].) The full picture of how the packet behaves is therefore 
quite complicated, though qualitatively it is still the front portion of the packet 
that is transmitted. The results of a 2-D simulation are shown in Figure (HJ 
and show a number of interesting features. For example, it is the part of the 
wavepacket that 'spins' into the barrier that is predominantly responsible for 
the transmitted wavepacket. The significance of the spin in the barrier region 
was clear from equation 1)7. 4|l . which showed that the spin vector generates a 
transverse current in the barrier region. These transverse currents are clearly 
displayed in Figure The motion near the barrier is highly complex, with the 
appearance of current loops suggesting the formation of vortices. Similar effects 
have been described by Hirschfelder et al. [HHI in the context of the Schrodinger 
theory. The streamline plots again show a slowing down in the barrier, which 
offsets the fact that it is the front of the packet that crosses the barrier. 

8 Spin Measurements 

We now turn to a second application of the local observables approach to quan- 
tum theory, namely to determine what happens to a wavepacket when a spin 
measurement is made. The first attempts to answer this question were made 
by Dewdney et al. |2H1 , who used the Pauli equation for a particle with zero 
charge and an anomalous magnetic moment to provide a model for a spin- 1/2 
particle in a Stern-Gerlach apparatus. Written in the STA, the relevant equation 
is 

dt^icjs. = - ^lB^a3, (8.1) 

2m 
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Figure 6: A 2-Dimensional simulation. The simulation uses a 2-D wavepacket 

in the x-y plane of width Ak = 0.04:A~^ and energy 5eV. The packet is incident 
perpendicular to the barrier, and the spin vector lies in the +z direction. The 
barrier is at x = and has width 2.5 A and height lOeV. The top two plots 
show streamlines in the plane. The top left plot shows streamlines about the 
bifurcation line, illustrating that the left side of the packet, which 'spins' into 
the barrier, is preferentially transmitted. The top right plot shows streamlines 
for a set of points near the front of the packet luith the same x value. Again, it 
is the left side of the packet that is transmitted. The reflected trajectories show 
quite complex behaviour, and in both plots the effect of the transverse current 
in the barrier is clear. The t-dependence of the streamlines in the top right 
plot is shown in the bottom two plots. The left-hand plot shows the streamlines 
in t-x space. Since the streamlines were started from the same x position and 
different y positions, the streamlines start from the same point and then spread 
out as the individual lines evolve differently. Again, it is possible to infer that 
the streamlines slow down as they pass through the barrier. The right-hand plot 
shows the t-y evolution of the same streamlines. In all plots distance is measured 
in A and time in units of 10~^'*s. 



and the current employed by Dewdney et al. is 

J = - — V($icr3$^). (8.2) 



m 



Dcwdncy et al. paramctcrisc the Pauh spinor 4> in terms of a density and three 
'Euler angles'. In the STA, this paramctcrisation takes the transparent form 
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where the rotor term is precisely that needed to parameterise a rotation in terms 
of the Euler angles. With this parameterisation, it is a simple matter to show 
that the current becomes 

J ^ — (VV^ + COS6IV0). (8.4) 
2 m 

But, as was noted in Section EITl the current defined by equation H8.2(l is not 
consistent with that obtained from the Dirac theory through a non-relativistic 
reduction. In fact, the two currents differ by a term in the curl of the spin 
vector IHEHI- 

To obtain a fuller understanding of the spin measurement process, an analy- 
sis based on the Dirac theory is required. Such an analysis is presented here. As 
well as dealing with a well-defined current, basing the analysis in the Dirac the- 
ory is important if one intends to proceed to study correlated spin measurements 
performed over spacelike intervals (i.e. to model an EPR-type experiment). To 
study such systems it is surely essential that one employs relativistic equations 
so that causality and the structure of spacetime are correctly built in. 

8.1 A Relativistic Model of a Spin Measurement 

As is shown in Section 14.31 the modified Dirac equation for a neutral particle 
with an anomalous magnetic moment [i, is 

V'0icr3 — = m-'iplo- (8-5) 

This is the equation we use to study the effects of a spin measurement, and it 
is not hard to show that equation H8.5|l reduces to l|8.1|l in the non-relativistic 
limit. Following Dewdney et al. pr2j we model the effect of a spin measurement 
by applying an impulsive magnetic field gradient, 

F = Bz5{t)i<T3. (8.6) 

The other components of B are ignored, as we are only modelling the behaviour 
of the packet in the z-direction. Around t — equation (|8.5|l is approximated 

by 

dt'ipias = Apz(5(t)73?A73, (8.7) 

where 

Ap = nB. (8.8) 
To solve (18 .711 we decompose the initial spinor ^pQ into 

tp^ = ^('/'o - 73V'o73), = ^(V'o + 73V'o73)- (8.9) 

Equation H8.7|l now becomes, for -0^ 

dti^^ = ApzS{t)ijh(T3 (8.10) 

with the opposite sign for t/jK The solution is now straightforward, as the 
impulse just serves to insert a phase factor into each of "0^ and i/;^: 
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If we now suppose that the initial tp consists of a positive-energy plane-wave 

^0 = L(p)$e*'^^(P-^ - ■^*) (8.12) 
then, immediately after the shock, is given by 

The spatial dependence of ip is now appropriate to two different values of the 
3- momentum, and , where 

p^=p + Apa3, p^=p — Ap(73. (8-14) 

The boost term L{p) corresponds to a different momentum, however, so both 
positive and negative frequency waves are required for the future evolution. 
After the shock, the wavefunction therefore propagates as 

^ = V'|e-*'^3P^-^ + + V+e"*'^^^^'^ + Vie*'^^^'^'^ (8.15) 

where 

p^jo^E^+p\ (8.16) 

p'^'jo^E'^ -p\ (8.17) 

and 

# = (m2+pT2)i/2^ (8 18) 

Both p^ and are defined similarly. 

Each term in H8.15|l must separately satisfy the free-particle Dirac equation, 
so it follows that 

pT^I = mibljo, (8.19) 
-pVI = m-ipljo, (8.20) 
which are satisfied together with 

TP^=i;l+Tpl. (8.21) 

The same set of equations hold for -0^. Dropping the arrows, we find that 

V'+ = ^(P7o0 + mV;) (8.22) 

ip- = ^(^ToV" - rn-ip) (8.23) 

which hold for both ■0^ and ip-^ . 

The effect of the magnetic shock on a monochromatic wave is to split the 
wave into four components, each with a distinct momentum. The positive fre- 
quency waves are transmitted by the device and split into two waves, whereas 
the negative frequency states are reflected. The appearance of the antiparticle 
states must ultimately be attributed to pair production, and only becomes sig- 
nificant for large B-fields. We examine this effect after looking at more physical 
situations. 
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8.2 Wavepacket Simulations 

For computational simplicity we take the incident particle to be localised along 
the field direction only, with no momentum components transverse to the field. 
This reduces the dimensionality of the problem to one spatial coordinate and the 
time coordinate. This was the set-up considered by Dewdney et al. and is 
sufficient to demonstrate the salient features of the measurement process. The 
most obvious difference between this model and a real experiment where the 
electron is moving is that, in our model, all four packets have group velocities 
along the field direction. 

The initial packet is built up from plane-wave solutions of the form 

^^e"''^/2$e''^3(p^-i?t) (8 24) 

which are superposed numerically to form a Gaussian packet. After the impulse, 
the future evolution is found from equation H8.15(l and the behaviour of the spin 
vector and the streamlines can be found for various initial values of The 
results of these simulations are plotted on the next few pages. 

In Figures |7| and |H1 we plot the evolution of a packet whose initial spin 
vector points in the ai direction ($ = exp{— icr27r/4}). After the shock, the 
density splits neatly into two equal-sized packets, and the streamlines bifurcate 
at the origin. As with the tunnelling simulations, we see that disjoint quantum 
outcomes are entirely consistent with the causal wavepacket evolution defined 
by the Dirac equation. The plot of the spin vector SA70 shows that immediately 
after the shock the spins are disordered, but that after a little time they sort 
themselves into one of the two packets, with the spin vector pointing in the 
direction of motion of the deflected packet. These plots are in good qualitative 
agreement with those obtained by Dewdney et al. |H2j , who also found that the 
choice of which packet a streamline enters is determined by its starting position 
in the incident wavepacket. 

Figure^lshows the results of a similar simulation, but with the initial spinor 
now containing unequal amounts of spin-up and spin-down components. This 
time we observe an asymmetry in the wavepacket split, with more of the density 
travelling in the spin- up packet. It is a simple matter to compute the ratio of the 
sizes of the two packets, and to verify that the ratio agrees with the prediction 
of standard quantum theory. 

As a final, novel, illustration of our approach, we consider a strong shock 
applied to a packet which is already aligned in the spin-up direction. For a 
weak shock the entire packet is deflected but, if the shock is sufficiently strong 
that the antiparticle states have significant amplitude, we find that a second 
packet is created. The significant feature of Figure [TUI is that the antiparticle 
states are deflected in the opposite direction, despite the fact that their spin 
is still oriented in the +z direction. The antiparticle states thus behave as if 
they have a magnetic moment to mass ratio of opposite sign. A more complete 
understanding of this phenomenon requires a field-theoretic treatment. The 
appearance of antiparticle states would then be attributed to pair production, 
with the antiparticle states having the same magnetic moment, but the opposite 
spin. (One of the crucial effects of the field quantisation of fermionic systems is 
to flip the signs of the charges and spins of antiparticle states.) 

The conclusions reached in this section are in broad agreement with those of 
Dewdney et ai. From the viewpoint of the local observables of the Dirac wave- 
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Figure 7: Splitting of a wavepacket caused by an impulsive B -field. The initial 
packet has a width of 1 x 10~^'^kg m in momentum space, and receives an 
impulse ofAp = 1 x 10~'^^kgms~^ . The top figure shows the probability density 
Jo at t = 0,1.3,2.6,3.9 x 10~^^s, with t increasing up the figure. The bottom 
figure shows streamlines in the (t, z ) plane. 

function (the current and spin densities), a Stern-Gerlach apparatus does not 

fulfil the role; of a classical measuring device. Instead, it behaves much more like 
a polariser, where the ratio of particles polarised up and down is dependent on 
the initial wavefunction. The B-field dramatically alters the wavefunction and 
its obscrvables, though in a causal manner that is entirely consistent with the 
predictions of standard quantum theory. The implications of these observations 
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Figure 8: Evolution of the relative spin-vector s A jo, in projection in the (x^z) 
plane. Immediately after the shock the spin-vectors point in all directions, but 
after about 2 x 10~"'^^s they sort themselves into the two packets, pointing in the 
+z and —z directions. 



for the interpretation of quantum mechanics are profound, though they are only 
slowly being absorbed by the wider physics community. (Some of these issues 
are debated in the collection of essays entitled ^Quantum Implications^ and 
in the recent book by Holland 

9 The Multiparticle STA 

So far we have dealt with the apphcation of the STA to single-particle quan- 
tum theory. In this section we turn to multiparticle theory. The aim here is to 
develop the STA approach so that it is capable of encoding multiparticle wave- 
functions, and describing the correlations between them. Given the advances in 
clarity and insight that the STA brings to single-particle quantum mechanics, 
we expect similar advances in the multiparticle case. This is indeed what we 
have found, although the field is relatively unexplored as yet. Here we highlight 
some areas where the multiparticle STA promises a new conceptual approach. 
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Figure 9: Splitting of a wavepacket with unequal mixtures of spin-up and spin- 
down components. The initial packet has $ = 1.618 — ia2, so more of the 
streamlines are deflected upwards, and the bifurcation point lies below the z = 
plane. The evolution of the spin vector is shown in the bottom plot. 



rather than attempting to reproduce the calculational techniques employed in 
standard approaches to many-body or many-electron theory. In particular, we 
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Figure 10: Creation of antiparticle states by a strong magnetic shock. The im- 
pulse used is Ap = 1 x 10^^^ kg ms^^ and the initial packet is entirely spin-up 
= 1). The packet travelling to the left consists of negative energy (antiparti- 
cle) states. 

concentrate on the unique geometric insights that the niultiparticle STA pro- 
vides — insights that are lost in the matrix theory. A prehminary introduction 
to the ideas developed here was given in 0J, though this is the first occasion 
that a full relativistic treatment has been presented. 

The 71-particle STA is created simply by taking n sets of basis vectors {7^}, 
where the superscript labels the particle space, and imposing the geometric 
algebra relations 



The fact that the basis vectors from distinct particle spaces anticommute means 
that we have constructed a basis for the geometric algebra of a 4n-dimensional 
configuration space. There is nothing uniquely quantum-mechanical in this 
idea — a system of three classical particles could be described by a set of 
three trajectories in a single space, or one path in a nine-dimensional space. 
The extra dimensions serve simply to label the properties of each individual 
particle, and should not be thought of as existing in anything other than a 
mathematical sense. This construction enables us, for example, to define a 
rotor which rotates one particle whilst leaving all the others fixed. The unique 
feature of the multiparticle STA is that it implies a separate copy of the time 
dimension for each particle, as well as the three spatial dimensions. To our 
knowledge, this is the first attempt to construct a solid conceptual framework 
for a multi-time approach to quantum theory. Clearly, if successful, such an 
approach will shed light on issues of locality and causality in quantum theory. 

The {7^} serve to generate a geometric algebra of enormously rich structure. 
Here we illustrate just a few of the more immediate features of this algebra. 



-yj^ii + liii = 0, i^ j 

lili + lll'i = 2?7^^ i = j. 



(9.1) 



These relations are summarised in the single formula 



(9.2) 



63 



It is our belief that the multiparticle STA will prove rich enough to encode 
all aspects of multiparticle quantum field theory, including the algebra of the 
fermionic creation/annihilation operators. 

Throughout, Roman superscripts are employed to label the particle space in 
which the object appears. So, for example, ip'^ and xjj'^ refer to two copies of the 
same 1-particle object ip, and not to separate, independent objects. Separate 
objects are given distinct symbols, or subscripts if they represent a quantity 
such as the current or spin-vector, which are vectors in configuration space 
with different projections into the separate copies of the STA. The absence of 
superscripts denotes that all objects have been collapsed into a single copy of 
the STA. As always, Roman and Greek subscripts are also used as frame indices, 
though this does not interfere with the occasional use of subscripts to determine 
separate projections. 



9.1 2-Particle Pauli States and the Quantum Correlator 

As an introduction to the properties of the multiparticle STA, we first consider 
the 2-particle Pauli algebra and the spin states of pairs of spin- 1/2 particles. 
As in the single-particle case, the 2-particle Pauli algebra is just a subset of the 
full 2-particle STA. A set of basis vectors is defined by 

= iho (9-3) 
^?=7ho^ (9.4) 

which satisfy 



12 1122 1221 2211 21 /nr^ 

f^i f^i = 7i 7o7i7o = 7i 7j7o7o = 7j7o7i 7o = ('j'^i ■ (9-5) 



So, in constructing multiparticle Pauli states, the basis vectors from differ- 
ent particle spaces commute rather than anticommute. Using the elements 
{Ijicrl, ial,iaj ia^} as a basis, we can construct 2-particle states. Here we have 
introduced the abbreviation 

ial = i'al (9.6) 

since, in most expressions, it is obvious which particle label should be attached 

to the i. In cases where there is potential for confusion, the particle label is put 
back on the i. The basis set {1, ia^, ia^, iaj ia\} spans a 16-dimensional space, 
which is twice the dimension of the direct product space of two 2-component 
complex spinors. For example, the outer-product space of two spin-1/2 states 
can be built from complex superpositions of the set 

1\ „ fl\ fO\ fl\ fl\ fo\ fo\ fo 



^ \0 



which forms a 4-dimensional complex space (8 real dimensions). The dimen- 
sionality has doubled because we have not yet taken the complex structure of 
the spinors into account. While the role of j is played in the two single-particle 
spaces by right multiplication by zcr| and ia^ respectively, standard quantum 
mechanics does not distinguish between these operations. A projection opera- 
tor must therefore be included to ensure that right multiplication by ia^ or ia^ 
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reduces to the same operation. If a 2-particle spin state is represented by the 
multivector 'ip, then ?/) must satisfy 

iPial = ipial (9.8) 

from which we find that 



(9.9) 



On defining 

E=l{l-iali4), (9.10) 

we find that 

E'^ = E (9.11) 

so right multiplication by is a projection operation. (The relation E'^ = E 
means that E is technically referred to as an 'idempotent' element.) It follows 
that the 2-particle state tp must contain a factor of E on its right-hand side. 
We can further define 

J = Eial = Eial = \{ial + iaf) (9.12) 

so that 

,f = -E. (9.13) 

Right-sided multiplication by J takes on the role of j for multiparticle states. 

The STA representation of a direct-product 2-particle Pauli spinor is now 
given by ip^cp'^E, where tp'^ and (p'^ are spinors (even multivectors) in their own 
spaces. A complete basis for 2-particle spin states is provided by 



E 

-ialE 
-ialE 
ial ia^E. 



(9.14) 



This procedure extends simply to higher multiplicities. All that is required 
is to find the 'quantum correlator' En satisfying 

Enia^s = E„ia^ = J„ for all j, k. (9.15) 

En can be constructed by picking out the j = 1 space, say, and correlating all 
the other spaces to this, so that 

n 

En = l[^{l-ialiai). (9.16) 

The value of En is independent of which of the n spaces is singled out and 
correlated to. The complex structure is defined by 

Jn = Eniai, (9.17) 
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where ia^^ can be chosen from any of the n spaces. To iUustrate this consider 
the case of n = 3, where 

E3 = \il-taltal){l^taltal) (9.18) 

= i(l - ial - i'^l - i^l (9-19) 

and 

J3 = 1(^0-3 + ial + *f^3 ^ i^l i^l i^l)- (9.20) 

Both i?3 and J3 are symmetric under permutations of their indices. 

A significant feature of this approach is that all the operations defined for 
the single-particle STA extend naturally to the multiparticle algebra. The rever- 
sion operation, for example, still has precisely the same definition — it simply 
reverses the order of vectors in any given multivector. The spinor inner prod- 
uct H3.12|l also generalises immediately, to 

(^-,0)5 = (i;„)"M(^^0> - (^^0Jn)ifT3], (9.21) 

where the right-hand side is projected onto a single copy of the STA. The factor 
of (En)^^ is included so that the state '1' always has unit norm, which matches 
with the inner product used in the matrix formulation. 



The Non-Relativistic Singlet State 

As an application of the formalism outlined above, consider the 2-particle singlet 
state |e), defined by 




(9.22) 



This is represented in the 2-particle STA by the multivector 

e = 7l(*^2 - ^^2)5(1 - *^3 ^^D- (9.23) 
The properties of e are more easily seen by writing 

e = i(l + ial *cr^)^(l + icTs 'i'^Dy^ial, (9.24) 

which shows how e contains the commuting idcmpotents ^(1 -f ia\ia2) and 
^(1 -|- ial ^'^i)- The normalisation ensures that 

^A{\{l + ialial)\{l + ialial)) 

= 1. (9.25) 

The identification of the idcmpotents in e leads immediately to the results 
that 

iale = \{ial - ial)\{l + ial ial)^ia\ = -ia^e (9.26) 

and 

iale = -ia^e, (9.27) 
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and hence that 

itrje = ial icr\e = — *CT2 — ^'^^^ ~ — icrje. (9.28) 

If M is an arbitrary even element in the Pauh algebra {M = M° + M^ial), it 
follows that e satisfies 

M^e = M^^e. (9.29) 

This now provides a novel demonstration of the rotational invariance of e. Under 
a joint rotation in 2-particle space, a spinor if) transforms to R^R^'ip, where 
and R^ are copies of the same rotor but acting in the two different spaces. From 
equation H9.29|l it follows that, under such a rotation, e transforms as 

e^-^ R^R'^e^ R^R^\^ e, (9.30) 
so that e is a genuine 2-particle scalar. 



Non-Relativistic Multiparticle Observables 

Multiparticle observables are formed in the same way as for single-particle states. 
Some combination of elements from the fixed {c^} frames is sandwiched between 
a multiparticle wavefunction ip and its spatial reverse V'^ ■ An important example 
of this construction is provided by the multiparticle spin-vector. In the matrix 
formulation, the fcth component of the particle- 1 spin vector is given by 

= (^la^l^) (9.31) 

which has the STA equivalent 

= -2"-i(jai).(VJV'^). (9.32) 

Clearly, the essential quantity is the bivector part of ij^J^^ which neatly gener- 
alises the single-particle formula. If we denote the result of projecting out from 
a multivector M the components contained entirely in the ith-particle space by 
{My, we can then write 

5^ = 2"-i(V^JV^)^ (9.33) 

The various subscripts and superscripts deserve some explanation. On both 
sides of equation (|9.33|) the superscript a labels the copy of the STA of inter- 
est. The subscript on the right-hand side as usual labels the fact that we are 
projecting out the grade-2 components of some multivector. The subscript a on 
the left-hand side is necessary to distinguish the separate projections of ipjip^. 
Had we not included the subscript, then and would refer to two copies 
of the same bivector, whereas S\ and S2 are different bivectors with different 
components. The reason for including both the subscript and the superscript 
on 5'" is that we often want to copy the individual bivectors from one space to 
another, without changing the components. 

We can hold all of the individual S"^ bivectors in a single multiparticle bivec- 
tor defined by 

5 = 2"-i(V'JV'^)2. (9.34) 
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Pauli 


Multivector 


Spin 


State 


Form 


Current 


1 TT) 


E2 


ia\ + «cr| 


ITi) 


-icTlE2 


ia\ — ial 


lit) 


—ia\E2 


—ial + *'^3 


1 ii) 




■ 1 ■ 9 

— ZtJg — ICTg 



Table 7: Spin Currents for 2-Particle Pauli States 

Under a joint rotation in n-particle space, tp transforms to Ri . . . Rni^ ^-nd S 
therefore transforms to 

i?i . . . R"SR''^ . . . = R^SlR^^ + ■■■ + R^'S^R''''. (9.35) 

Each of the separate projections of the spin current is therefore rotated by the 
same amount, in its own space. That the definition (|9.34l) is sensible can be 
checked with the four basis states (|9.14() . The form of S for each of these 
is contained in Table Multiparticle observables for the 2-particle case are 
discussed further below. 

Other observables can be formed using different fixed multivectors. For ex- 
ample, a 2-particle invariant is generated by sandwiching a constant multivector 
S between the singlet state e, 

M = eY,e^. (9.36) 

Taking S = 1 yields 

M = ee^ = 2i(l + ia^ ial)\{l + ial '^'^l) 

= ^(1 -I- ia\ ia\ -t- ia\ icl -f 4173 icr|), (9.37) 

which rearranges to give 

ialial = 2ee'^ -\. (9.38) 

This equation contains the essence of the matrix result 

at,, al,,= 251,5",, -61,51, (9.39) 

where a, b, a', h' label the matrix components. This matrix equation is now 
seen to express a relationship between 2-particle invariants. Further invariants 
are obtained by taking S = i^i^, yielding 

M = ei^i^e^ = \ {i^i^ + a\ af + a^ aj + al a^). (9.40) 

This shows that both iaj. iaf, and a], a'^ are invariants under 2-particle rotations. 
In standard quantum mechanics these invariants would be thought of as arising 
from the "inner product" of the spin vectors a} and af. Here, we have seen 
that the invariants arise in a completely different way by looking at the full 
multivector ee^. 

The contents of this section should have demonstrated that the multiparticle 
STA approach is capable of reproducing most (if not all) of standard multipar- 
ticle quantum mechanics. One important result that follows is that the unit 
scalar imaginary j can be completely eliminated from quantum mechanics and 
replaced by geometrically meaningful quantities. This should have significant 
implications for the interpretation of quantum mechanics. 
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9.2 Comparison with the 'Causal' Approach to 
Non-Relativistic Spin States 

As an application of the techniques outhned above, we look at the work of 
Holland on the 'Causal interpretation of a system of two spin-1/2 particles' 
This work attempts to give a non-relativistic definition of local observables in 
the higher-dimensional space of a 2-particle wavefunction. As we have seen, such 
a construction appears naturally in our approach. Holland's main application 
is to a Bell inequality type experiment, with spin measurements carried out on 
a system of two correlated spin-1/2 particles by Stern-Gerlach experiments at 
spatially separated positions. Such an analysis, though interesting, will only be 
convincing if carried out in the fully relativistic domain, where issues of causality 
and superluminal propagation can be coherently addressed. We intend to carry 
out such an analysis in the future, using the STA multiparticle methods, and 
the work below on the observables of a 2-particle system can be seen as part of 
this aim. 

The aspect of Holland's work that concerns us here (his Section 3 and Ap- 
pendix A) deals with the joint spin-space of a system of two non-relativistic 
spin-1/2 particles. The aim is to show that 'all 8 real degrees of freedom in 
the two body spinor wavefunction may be interpreted (up to a sign) in terms 
of the properties of algebraically interconnected Euclidean tensors' [Hi]- Hol- 
land's working is complex and requires a number of index manipulations and 
algebraic identities. Furthermore, the meaning of the expressions derived is far 
from transparent. Using the above techniques, however, the significant results 
can be derived more efficiently and in such a way that their geometric meaning 
is made much clearer. Rather than give a line-by-line translation of Holland's 
work, we simply state the key results in our notation and prove them. 

Let V' = ipE be a 2-particle spinor in the correlated product space of the one- 
particle spin-spaces (the even subalgebra of the Pauh algebra). The observables 
of this 2-particle system are formed from projections of bilinear products of the 
form il^Tij}, where F is an element of the 2-particle Pauli spinor algebra. For 
example, the two 3-dimensional spin- vectors, s}, 83, are defined by 

isl+isj = 2^Jil) (9.41) 

where the right-hand side can be written in the equivalent forms 

ipj^jj = Tjjialijj = i^iali). (9.42) 

The formula 1)9. 41|) is a special case of equation H9.34() where, as we are working 
in a 2-particle system, the projection onto bivector parts is not required. The 
vectors s\ and correspond to the two spin vectors defined by Holland, with 
the explicit correspondence to his Sik and S2k given by 

S^k = -i^sl)-{ial), S2k^-{^sl)-{ial). (9.43) 

An important relation proved by Holland is that the vectors s\ and are 
of equal magnitude, 

{sir = isjr = 2n-p^, (9.44) 

where 

p = 2{i/jiP^) (9.45) 
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and an explicit form for is to be determined. To prove this result, we write 
the formulae for the components of sj and (19.431) in the equivalent forms 

Sik = -2{i;Ual,p) ■ (tal), S^k - -2(^^(7^ V-) • (9.46) 

But, in both cases, the term tpHa'^tp contains a bivector sandwiched between 
two idempotents, so is of the form E . . . E. This sandwiching projects out the 
ial and ia^ components of the full bivector, and ensures that that the terms 
have equal magnitude. The inner products in (|9.46f) can therefore be dropped 
and we are left with 

^Uat^ = SakJ (9.47) 

where the a = 1,2 labels the two separate spin- vectors. It follows immediately 
from equation H9.47|l that 

(S^)2 = -2(^t,^a^^t,^a^)^ (9.48) 

where the a's are not summed. But the quantity V'V'^ contains only scalar and 
4-vector components, so we find that 

iali^tp'' ial = iali^ipUal = -3(?AV'^) + (V'V'^)4 = V"/-^ - 2p. (9.49) 

(This result follows immediately from Ukiaak = — m, which is valid for any 
vector a in the single-particle Pauli algebra.) Inserting equation (|9.49|) back 
into (|9.48f) we can now write 

= (s2)'-2p2_2(V'^t^^t)^ (9.50) 

which shows that s\ and s\ are indeed of equal magnitude, and enables us to 
identify f2 as 

i(3p2-2(V'V'^VV'^)). (9.51) 

In addition to p, s\ and Holland defines a tensor Sjk whose components 
are given by 

S,k = -2{'il;ha]i<7l'il;). (9.52) 
This object has the simple frame-free form 

T = (V'V^)4. (9.53) 

Between them, p, sj, s\ and T pick up 7 of the possible 8 degrees of freedom 
in ij). The remaining freedom lies in the phase, since all of the observables 
defined above are phase- invariant. Encoding this information caused Holland 
some difficulty, but in the STA the answer is straightforward, and is actually 
already contained in the above working. The crucial observation is that as 
well as containing only scalar and 4-vector terms, the quantity ^/^^-^ is invariant 
under rotations. So, in addition to the scalar p, the 4-vector components of '0^V' 
must pick up important rotationally-invariant information. Furthermore, since 
is of the form E . . . E, the four-vector component of V'^V' contains only two 
independent terms, which can be taken as a complex combination of the i(T2 ?cr| 
term. This is seen most clearly using an explicit realisation of ijj. Suppose that 
we write 

= [p — i<7\q — i<y2^' + i<y\ ia2s)E (9.54) 
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where p, q, r, and s are complex combinations of 1 and J, we then find that 

V^t^ = [p + 2ial ialips ~ qr)]E. (9.55) 

This shows expHcitly that the additional complex invariant is given by ps — qr. 
It is this term that picks up the phase of "0, and we therefore define the complex 
quantity 

a — {ipicr2 i<J2ilJ^) — {ipi(j\ zcr|?A^)icr3, (9.56) 

which is the STA equivalent of the complex scalar p defined by Holland. The 
complex scalar a is invariant under rotations, and under the phase change 

^ ije'^'t' (9.57) 

a transforms as 

at^ae^'^*^^ (9.58) 

The set {p, a, s}, T} encode all the information contained in the 2-particle 
spinor ij:, up to an overall sign. They reproduce the quantities defined by Hol- 
land, but their STA derivation makes their properties and geometric origin much 
clearer. 



9.3 Relativistic 2-Particle States 

The ideas developed for the multi-particle Pauli algebra extend immediately 
to the relativistic domain. The direct product of the two single-particle spinor 
spaces (the even subalgebras) now results in a space of 8 x 8 = 64 real dimen- 
sions. Unlike the single-particle case, this space is not equivalent to the even 
subalgebra of the full 8-dimensional algebra. The full algebra is 256-dimcnsional, 
and its even subalgebra is therefore 128-dimensional. It is not yet clear whether 
the remaining 64-dimensional space which is not picked up by sums of direct- 
product states could be of use in constructing 2-particle wavef unctions, and 
for the remainder of this section we work only with the space obtained from 
sums of direct-product states. Post-multiplying the direct-product space by the 
quantum correlator E reduces it to 32 real dimensions, which are equivalent to 
the 16 complex dimensions employed in standard 2-particle relativistic quan- 
tum theory. All the single-particle observables discussed in Section rOl extend 
simply. In particular, we define the vectors 

J ^{i'ill+ll)^)! (9.59) 
s = (^(73+7lWi (9-60) 

which are respectively the 2-particle current and spin-vector. (The calligraphic 
symbol J is used to avoid confusion with the correlated bivector J.) We also 
define the spin bivector S by 

S^{iPJi!)2. (9.61) 

Of particular interest are the new Lorentz-invariant quantities that arise in 
this approach. From the work of the preceding section, we form the quantity 
which decomposes into 

= (V^V')o,8 + (V^V')4. (9.62) 
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The grade-0 and grade-8 terms are the 2-particle generahsation of the scalar 
+ pseudoscalar combination — pexp(i/3) found at the single-particle level. 
Of greater interest are the 4-vector terms. These offer a wealth of Lorentz- 
invariant 2-particle observables, the meaning of which we are only beginning to 
appreciate. Such invariants are rarely seen in the traditional matrix approach. 



The Relativistic Singlet State and Relativistic Invariants 

Our task here is to find a relativistic analogue of the Pauli singlet state discussed 
in Section IHTI Recalling the definition of e ()9.23|) . the property that ensured 
that e was a singlet state was that 

iale^-iale, /c = 1...3. (9.63) 

In addition to H9.63|l a relativistic singlet state, which we will denote as 77, must 
satisfy 

'Jlv = -alv, k = 1...3. (9.64) 

It follows that rj satisfies 

zSy = a\a\al-q = -cricrlcrlv = ^^V (9-65) 

so that 

77 = —i^f'rj (9.66) 

=^ ri=\{l-i^i^)ri. (9.67) 

The state 77 can therefore be constructed by multiplying e by the idempotent 
^(1 — i^i"^). We therefore define 

which is normalised such that (77,77)5 — 1. The invariant 77 satisfies 

ialr] = iale\{\-i^i^) = -iolri A: = 1 . . . 3 (9.69) 

and 

(^Iri^-^^li^i'^'n^i'^icFl-n^-c^l-n fc = 1...3. (9.70) 
These results are summarised by 

MI77 = M'^'q, (9.71) 

where M is an even multivector in either the particle-1 or particle-2 STA. The 
proof that 7; is a relativistic invariant now reduces to the simple identity 

R^R^i] = R^R^r] = 77, (9.72) 

where i? is a single-particle relativistic rotor. 

Equation H9.71|l can be seen as arising from a more primitive relation between 
vectors in the separate spaces. Using the result that 7o7o commutes with rj, we 
can derive 

ilvio = ilililvihhl 
222 2 

= 7'^7o, (9.73) 
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and hence wc find that, for an arbitrary vector a, 

a'vio = a^ml (9.74) 
Equation H9.71|l now foUows immediately from H9.74|l by writing 

{abfr, = ah'r^^Wo 

7,2 2 2 2 

= [hafr^. (9.75) 

Equation (|9.74|l can therefore be viewed as the fundamental property of the 
relativistic invariant rj. 

From rj a number of Lorentz-invariant 2-particle multivectors can be con- 
structed by sandwiching arbitrary multivectors between and i). The simplest 
such object is 

r/f) — e^{l — i^i'^)e 

i(l + icrl ial + ia\ iaj + icrl icr^^il - i^i'^) 
^^{l-^^^^)-^{alal-^al^al). (9.76) 

This contains a scalar + pseudoscalar (grade-8) term, which is obviously invari- 
ant, together with the invariant grade-4 multivector (cr^ cr^ — icr^ The next 
simplest object is 

= i(7o7o-7hfe)(l-^'^')- (9.77) 
On defining the bivector 

K = 717^2 (9 

and the 2-particle pseudoscalar 

W = i^i^=i^i^ (9.79) 

the invariants from H9.77|l are simply K and WK. That W is invariant under 
rotations is obvious, and the invariance of K under joint rotations in the two 
particle spaces follows from equation (|9.77|l . The bivector K is of the form of a 
'doubling' bivector discussed in [^2], where such bivectors are shown to play an 
important role in the bivector realisation of many Lie algebras. 
From the definition of K H9.78|l . wc find that 

KAK = ~2^hlihl + (7hfe)A(7l7') 

= 2{alal~iaUal), (9.80) 

which recovers the grade-4 invariant from H9.76|l . The full set of 2-particle 
invariants constructed from K are summarised in Table [S] These invariants 
are well-known and have been used in constructing phenomenological models of 
interacting particles |f)3l Ifi4| . The STA derivation of the invariants is quite new, 
however, and the fundamental role played by the bivector K is hidden in the 
matrix formalism. 
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Type of 




Invariant 


Interaction 


Grade 


1 


Scalar 





K 


Vector 


2 


KAK 


Bivector 


4 


WK 


Pseudovector 


6 


W 


Pseudoscalar 


8 



Table 8: 2- Particle Relativistic Invariants 



9.4 Multiparticle Wave Equations 

In order to extend the local-observables approach to quantum theory to the 
multiparticle domain, we need to construct a relativistic wave equation satisfied 
by an n-particle wavefunction. This is a subject that is given little attention 
in the literature, with most textbooks dealing solely with the field-quantised 
description of an n-particle system. An n-particle wave equation is essential, 
however, if one aims to give a relativistic description of a bound system (where 
field quantisation and perturbation theory on their own are insufficient). A 
description of this approach is given in Chapter 10 of Itzykson & Zuber |37j . 
who deal mainly with the Bethe-Salpeter equation for a relativistic 2-particle 
system. Written in the STA, this equation becomes 

(jyi - mi)(jV2 - m2Mr, s) = I{r, s)V'(r, s) (9.81) 

where j represents right-sided multiplication by J, /(r, s) is an integral operator 
representing the inter-particlc interaction, and 

V'-T^T^, ^'^lll^, (9.82) 

with r and s the 4-D positions of the two particles. Strictly, we should have 
written V,^. and instead of simply and V^. In this case, however, the 
subscripts can safely be ignored. 

The problem with equation H9.81|l is that it is not first-order in the 8- 
dimensional vector derivative V = -I- V^. We are therefore unable to gen- 
eralise many of the simple first-order propagation techniques discussed in Sec- 
tion Clearly, we would like to find an alternative to (|9.81(l which retains 
the first-order nature of the single-particle Dirac equation. Here we will simply 
assert what we believe to be a good candidate for such an equation, and then 
work out its consequences. The equation we shall study, for two free spin-1/2 
particles of masses mi and m2 respectively, is 

[Zi + Yl\ ^(a;) (j^i + ^^2) = 2^(x). (9.83) 

We can assume, a priori, that ip is not in the correlated subspace of the the 
direct-product space. But, since E commutes with 173 + 17!, any solution 
to (|9.83(l can be reduced to a solution in the correlated space simply by right- 
multiplying by E. Written out explicitly, the vector x in equation (|9.83(l is 

x = r^ ^ -flrf" + jlsf" (9.84) 
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where {r'^,s'^} are a set of 8 independent components for ip. Of course, all 
particle motions ultimately occur in a single space, in which the vectors r and 
s label two independent position vectors. We stress that in this approach there 
are two time-like coordinates, r*^ and s°, which is necessary if our 2-particle 
equation is to be Lorentz covariant. The derivatives and are as defined 
by equation H9.82|) , and the 8-dimensional vector derivative V = is given by 

V = Vi+V^ (9.85) 

Equation ipH?^ can be derived from a Lorentz-invariant action integral in 8- 
dimensional configuration space in which the f /mi and 1 / m2 factors enter via 
a linear distortion of the vector derivative V. We write this as 

— + — U/i(V), (9.86) 
mi m2 / 

where h is the linear mapping of vectors to vectors defined by 

hia)^l-(a-r')ll + ^^{a-r')ll. (9.87) 

This distortion is of the type used in the gauge theory approach to gravity 
developed in Jl^l 1651 1661 l?)7j , and it is extremely suggestive that mass enters 
equation H9.83|l via this route. 

Any candidate 2-particle wave equation must be satisfied by factored states 
of the form 

^ = <t'Hr')xHs')E, (9.88) 
where (j)^ and a-re solutions of the separate single-particle Dirac equations, 

V0 = -mi(/)i73, Vx = -TO2xi73- (9.89) 

To verify that our equation (|9.83|l meets this requirement, we substitute in the 
direct-product state (|9.88() and use H9.89|l to obtain 

fYl + Z!") ^^^^E (77! + Z7|) = ~<ly\'E (171 + tji) (171 + 17!) , (9.90) 
\mi m2 / 

where we have used the result that commutes with 0^. Now, since 773 and 
77! anticommute, we have 

H + nl) H + - -2 (9.91) 

so that ^ ^ ^ 

— + — ) <I>\^E {ijl + ijI) = 2c^'x^E, (9.92) 

and H9.83|l is satisfied. Equation (|9.83() is only satisfied by direct-product states 
as a result of the fact that vectors from separate particle spaces anticommute. 
Hence equation H9.83(l does not have an equivalent expression in terms of the 
direct-product matrix formulation, which can only form commuting operators 
from different spaces. 
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9.5 The Pauli Principle 

In quantum theory, indistinguishable particles must obey either Fermi-Dirac or 

Bosc-Einstcin statistics. For fcrmions this requirement results in the Pauli exclu- 
sion principle that no two particles can occupy a state in which their properties 
are identical. At the relativistic multiparticle level, the Pauli principle is usu- 
ally encoded in the anticommutation of the creation and annihilation operators 
of fermionic field theory. Here we show that the principle can be successfully 
encoded in a simple geometrical manner at the level of the relativistic wavefunc- 
tion, without requiring the apparatus of quantum field theory. 
We start by introducing the gradc-4 multivector 

/ = roFirsrs, (9.93) 

where 

^^ = ^U + ^l)- (9-94) 
It is a simple matter to verify that / has the properties 

= -1, (9.95) 

and 

Hi = -rl = ^l- (9-96) 

It follows that / functions as a geometrical version of the particle exchange 
operator. In particular, acting on the 8-dimensional position vector x = + 
we find that 

Ixl = r^ + (9.97) 

where 

r2 = = jy. (9.98) 

So / can certainly be used to interchange the coordinates of particles 1 and 2. 
But, if I is to play a fundamental role in our version of the Pauli principle, we 
must first confirm that it is independent of our choice of initial frame. To see 
that it is, suppose that we start with a rotated frame {i?7^.R} and define 

= i= [r\R' + RW^ie) = R^R^T,R?R\ (9.99) 

The new F^ give rise to the rotated 4- vector 

I' = R^R^IR^RK (9.100) 

But, acting on a bivector in particle space 1, we find that 

la^Ab^I = -{Ia^I)A{Ib^I) = -a^Ab^, (9.101) 

and the same is true of an arbitrary even element in either space. More generally, 
I . . .1 applied to an even element in one particle space flips it to the other 
particle space and changes sign, while applied to an odd element it just flips the 
particle space. It follows that 

IR^R^ = R^IR^ = R^R^I, (9.102) 
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and substituting this into H9.100|l wc find that /' = /, so / is indeed independent 
of the chosen orthonormal frame. 

We can now use the 4-vector / to encode the Pauh exchange principle ge- 
ometrically. Let ip{x) be a wavefunction for two electrons. Our suggested rel- 
ativistic generalization of the Pauli principle is that i'ix) should be invariant 
under the operation 

V'(x) ^ Ii>{IxI)I. (9.103) 

For n-particle systems the extension is straightforward: the wavefunction must 
be invariant under the interchange enforced by the /'s constructed from each 
pair of particles. 

We must first check that (|9.103|) is an allowed symmetry of the 2-particle 
Dirac equation. With x' defined as Ixl it is simple to verify that 

V^r' = v.? + V] = /V/, (9.104) 

and hence that 

V = /V^'/. (9.105) 

So, assuming that satisfies the 2-particle equation (|9.83(l with equal masses 
m, we find that 

V[mixI)I] {i^l + Z7|) = -IV^Mx'V (*73 + nl) (9.106) 
= mlil}{x') {i-il -f ^7^) / (^73^ -t- i-il) . 

But «73 -|- «7| is odd and symmetric under interchange of its particle labels. It 
follows that 

^ (^73 + nl) ^ = (*73 + nl) (9.107) 

and hence that 

V[Ii;{IxI)I] {ij^ + i-fi) 2mIi}j{IxI)I. (9.108) 

So, if tl){x) is a solution of the 2-particle equal-mass Dirac equation, then so to 
is Ii>{IxI)I. 

Next we must check that the proposed relativistic Pauli principle deals cor- 
rectly with well-known elementary cases. Suppose that two electrons are in the 
same spatial state. Then we should expect our principle to enforce the condition 
that they are in an antisymmetric spin state. For example, consider ia\ — icr|, 
the spin singlet state. We find that 

I{ial-i(jl)I = -ial + ial, (9.109) 

recovering the original state, which is therefore compatible with our principle. 
On the other hand 

I{ial+ial)I = -{ial+ial), (9.110) 

so no part of this state can be added in to the wavefunction, which again is 
correct. In conclusion, given some 2-particle solution ili{x), the corresponding 
state 

^ ipix) + I^{IxI)I (9.111) 

still satisfies the Dirac equation and is invariant under ^pi^) ^ I'ip{IxI)I. We 
therefore claim that the state tpi is the correct relativistic generalisation of a 
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state satisfying the Pauli principle. In deference to standard quantum theory, 
we refer to equation IjO. 111(1 as an antisymmetrisation procedure. 

The final issue to address is the Lorentz covariance of the antisymmetrisation 
procedure 1(9.111(1 . Suppose that we start with an arbitrary wavefunction 'ipi^) 
satisfying the 2-particle equal- mass equation ((9.83(1 . If we boost this state via 

^{x) ^ ^'{x) = R^R^^PiR^R^xR^R^) (9.112) 

then ip'{x) also satisfies the same equation ((9.83(1 . The boosted wavefunction 
ip'{x) can be thought of as corresponding to a different observer in relative 
motion. The boosted state ip'i^) can also be antisymmetrised to yield a solution 
satisfying our relativistic Pauli principle. But, for this procedure to be covariant, 
the same state must be obtained if we first antisymmetrise the original ^'{x), 
and then boost the result. Thus we require that 

Sti^iSxS) + isij{isxsi)i = s[i^{SxS) + ii;{sixis)i] (9.113) 

where S = R^R^. Equation ((9.113(1 reduces to the requirement that 

IStpilSxSI) = Sli^iSIxIS) (9.114) 
which is satisfied provided that 

IS = SI (9.115) 

or 

R^R^I = IR^R^. (9.116) 

But we proved precisely this equation in demonstrating the frame-invariance of 
/, so our relativistic version of the Pauli principle is Lorentz invariant. This 
is important as, rather like the inclusion of the quantum correlator, the Pauli 
procedure discussed here looks highly non-local in character. 

9.6 8-Dimensional Streamlines and Pauli Exclusion 

For a single Dirac particle, a characteristic feature of the STA approach is that 
the probability current is a rotated/dilated version of the 70 vector, J = ipjo'ip- 
This current has zero divergence and can therefore be used to define streamlines, 
as discussed in Section|Z| Here we demonstrate how the same idea extends to the 
2-particle case. We find that the conserved current is now formed from -0 acting 
on the 7q -I-70 vector, and therefore exists in 8-dimensional configuration space. 
This current can be used to derive streamlines for two particles in correlated 
motion. This approach should ultimately enable us to gain a better insight into 
what happens in experiments of the Bell type, where spin measurements on 
pairs of particles are performed over spacelike separations. We saw in Section|Hl 
how the local observables viewpoint leads to a radical re-interpretation of what 
happens in a single spin-measurement, and we can expect an equally radical 
shift to occur in the analysis of spin measurements of correlated particles. As 
a preliminary step in this direction, here we construct the current for two free 
particles approaching each other head-on. The streamlines for this current are 
evaluated and used to study both the effects of the Pauli antisymmetrisation and 
the spin-dependence of the trajectories. This work generalises that of Dewdney 
et al. [681 12 7[ to the relativistic domain. 
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We start with the 2-particle Dirac equation H9.83|l . and multiply on the right 
by E to ensure the total wavefunction is in the correlated subspace. Also, since 
we want to work with the indistinguishable case, we assume that both masses 
are m. In this case our basic equation is 

V^IjE {ij^ + i"fi) = 2m^pE (9.117) 

and, since 

i?(«73+*73) =^(7o+7o), (9.118) 
equation H9.117fl can be written in the equivalent form 

VipE (70^ + 7^) = -2mV'J. (9.119) 
Now, assuming that satisfies ip = ipE, we obtain 

VV'(7o + lo)i' = -2mVj J^A, (9.120) 
and adding this equation to its reverse yields 

V^(7o + la)^ + ^(7o + 7o)'/'V = 0. (9.121) 
The scalar part of this equation gives 

V-(^(7o+7oWi = 0, (9.122) 
which shows that the current we seek is 

:^=(^(7o+7o)^)i, (9.123) 

as defined in equation 1)9.59(1 . The vector J' has components in both particle- 1 
and particle-2 spaces, which we write as 

J--Jl + Jl (9.124) 

The current is conserved in ei^/ii-dimensional space, so its streamlines never 
cross there. The streamlines of the individual particles, however, are obtained 
by integrating Ji and J2 in ordinary 4-d space, and these can of course cross. 
An example of this is illustrated in Figure lTTl which shows the streamlines corre- 
sponding to distinguishable particles in two Gaussian wavepackets approaching 
each other head-on. The wavefunction used to produce this figure is just 

^ = cj,\r^)x\s^)E, (9.125) 

with (j) and x being Gaussian wavepackets, moving in opposite directions. Since 
the distinguishable case is assumed, no Pauli antisymmetrisation is used. The 
individual currents for each particle are given by 

Ji(r,s) = (j){r)-ia(j){r) (x(s)x(s)), 

J2(r, s) = x(s)7ox(s) (0(r)0(r)) (9.126) 

and, as can be seen, the streamlines (and the wavepackets) simply pass straight 
through each other. 
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Distinguishable Particies 




Figure 11: Streamlines generated by the unsymmetrised 2-particle wavefunction 
ip = (f>^{r^)x'^is'^)E. Time is shown on the vertical axis. (f> and x o,re Gaussian 
wavepackets moving in opposite directions and the 'collision' is arranged to take 
place at t = 0. The lack of any antisymmetrisation applied to the wavefunction 
means that the streamlines pass straight through each other. 

An interesting feature emerges in the individual currents in H9.126I) . One of 
the main problems with single-particle Dirac theory is that the current is al- 
ways positive-definite so, if we wish to interpret it as a charge current, it fails to 
represent antiparticles correctly. The switch of sign of the current necessary to 
represent positrons is put into conventional theory essentially 'by hand', via the 
anticonimutation and normal ordering rules of fermionic field theory. In equa- 
tion (|9.126() . however, the norm (xx) of the second state multiplies the current 
for the first, and vice versa. Since (xx) can be negative, it is possible to obtain 
currents which fiow backwards in time. This suggests that the required switch of 
signs can be accomplished whilst remaining wholly within a wavefunction-based 
approach. An apparent problem is that, if only one particle has a negative norm 
state — say for example x has (xx) < — then it is the 4> current which is 
reversed, and not the x current. However, it is easy to see that this objection 
is not relevant to indistinguishable particles, and it is to these we now turn. 

We now apply the Pauli symmetrization procedure of the previous subsec- 
tion to the wavefunction of equation l|9.125l) . so as to obtain a wavefunction 
applicable to indistinguishable particles. This yields 



from which we form J7i and J72, as before. We must next decide which spin states 
to use for the two particles. We first take both particles to have their spin vectors 
pointing in the positive z-direction, with all motion in the iz-direction. The 
resulting streamlines are shown in Figure [T^ (top) . The streamlines now 'repel' 
one another, rather than being able to pass straight through. The corrugated 
appearance of the lines near the origin is the result of the streamlines having 
to pass through a region of highly-oscillatory destructive interference, since the 
probability of both particles occupying the same position (the origin) with the 
same spin state is zero. If instead the particles are put in different spin states 




(9.127) 
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Spins Aligned 




Spins Anti-aligned 




Figure 12: Streamlines generated by the antisymmetrised 2-particle wavefunction 
Tp — [(l)^{r^)x'^is'^) ~ X^{''''^)4''^is^)]E!. The individual wavepackets pass through 
each other, but the streamlines from separate particles do not cross. The top 
figure has both particles with spins aligned in the +z-direction, and the bottom 
figure shows particles with opposite spins, with (p in the +z direction, and x ^n 
the ~z direction. Both wavepackets have energy 527KeV and a spatial spread of 
~ 20pm. The spatial units are lO^^^m and the units of time are 10~^^s. The 
effects of the antisymmetrisation are only important where there is significant 
wavepacket overlap. 



then the streamUnes shown at the bottom of Figure El resuh. In this case 
there is no destructive interference near the origin, and the streamUnes are 
smooth there. However, they still repel! The explanation for this lies in the 
symmetry properties of the 2-particle current. Given that the wavefunction tp 
has been antisymmetrised according to our version of the Pauli principle, then 
it is straightforward to show that 

IJ{IxI)I = J{x). (9.128) 
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It follows that at the same spacetime position, encoded by Ixl = a; in the 2- 
particle algebra, the two currents J7i and J^2 are equal. Hence, if two streamlines 
ever met, they could never separate again. For the simulations presented here, 
it follows from the symmetry of the set-up that the spatial currents at the origin 
are both zero and so, as the particles approach the origin, they are forced to 
slow up. The delay means that they are then swept back in the direction they 
have just come from by the wavepacket travelling through from the other side. 
We therefore see that 'repulsion' as measured by streamlines has its origin in 
indistinguishability, and that the spin of the states exerts only a marginal effect. 



10 Further Applications 

In this section we briefly review two further applications of spacetime algebra to 
important areas of electron physics. The first of these is classical and semiclassi- 
cal mechanics. As well as simplifying many calculations in quantum mechanics, 
spacetime algebra is well suited to handling problems in classical mechanics 
where electrons are treated as point charges following a single trajectory. In 
recent years there has been considerable interest in finding modifications to the 
simple classical equations to include the effects of spin, without losing the idea 
of a definite trajectory |^|^. One of the aims behind this work is to find a 
suitable classical model which can be quantised via the path- integral route |70| . 
One of the more promising candidates is discussed here, and we outline some 
improvements that could repair some immediate defects. 

The second application discussed here is to Grassmann algebra and the asso- 
ciated 'calculus' introduced by Berezin |71) . Grassmann quantities are employed 
widely in quantum field theory, and the Berezin calculus plays a crucial role in 
the path-integral quantisation of fermionic systems. Here we outline how many 
of the calculations can be performed within geometric algebra, and draw atten- 
tion to some work in the literature. We do not attempt a more detailed analysis 
of path-integral quantisation here. 



10.1 Classical and Semiclassical Mechanics 

The Lorentz force law for a point-particle with velocity w, mass m and charge q 
is 

v = ^F-v (10.1) 
m 

where v denotes differentiation with respect to the affine parameter and F is 
the external electromagnetic field bivector. Any future-pointing unit timelike 
vector can be written in terms of a rotor acting on a fixed vector 70, 

v{t) = R{t)joR{t) (10.2) 

from which we find 

v = {2RR)-v. (10.3) 
The quantity RR is a bivector, so we can recover equation (jl0.1|l by setting 

2RR = —F (10.4) 

m 

=J> R= ^FR. (10.5) 

2m ^ ' 



82 



This is not the only equation for R that is consistent with 1)11). l|l . since any 
bivector that commutes with v could be added to F. However, (110. 5() is without 
doubt the simplest equation available. 

It turns out that equation 1)10. 5|l is often easier to analyse than (|10.1() . as 
was first shown by Hestenes jI4] . Furthermore, we can extend this approach to 
include a classical notion of spin. Let us suppose that, as well as describing the 
tangent vector u, the rotor R determines how a frame of vectors is transported 
along the curve. We can then define the spin vector as the unit spatial vector 

s = RjsR, (10.6) 

which matches the definition given for the quantum observable. If we now 
assume that equation 1)10. 5|l is valid, we find that the spin-vector satisfies the 
equation 

s^—F-s, (10.7) 
m 

which gives the correct precession equation for a particle of gyromagnetic ratio 
2 It follows that g = 2 can be viewed as the natural value from the 

viewpoint of the relativistic classical mechanics of a rotating frame — a striking 
fact that deserves to be more widely known. We can use the same approach to 
analyse the motion of a particle with a (^-factor other than 2 by replacing (|10.5() 
with 

R^^[FR+ig/2-l)RBi (10.8) 

which reproduces the Bargmann- Michel- Telegdi equation employed in the anal- 
ysis of spin- precession measurements 48^ . 

A remarkable aspect of the Dirac theory is that the current iplo^' the 
momentum (which is defined in terms of the momentum operator) are not neces- 
sarily colinear. This suggests that a more realistic classical model for an electron 
should employ an independent quantity for the momentum which is not neces- 
sarily related to the tangent vector to the spacetime trajectory. Such a model 
was proposed by Barut & Zanghi who did not employ the STA, and was 
analysed further in (see also [Z21)- Written in the STA, the action proposed 
by Barut & Zanghi takes the form 

S = J dX {■tpias'il) + p{i ~ ■tp-fo'tjj) + qA{x)'tp-foi') (10.9) 

where the dynamical variables are a;(A), p(A) and tpW- Variation with respect 
to these variables yields the equations 

X = i>-fa^ (10.10) 

P^qF-x (10.11) 

Tpias = P^lo (10.12) 

where 

P^p-qA. (10.13) 

These constitute a set of first-order equations so, with x, p and i}} given for 
some initial value of A, the future evolution is uniquely determined. Equations 
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(|10.1U|1 - (|10.12|I contain a number of unsatisfactory features. One does not ex- 
pect to see P entering the Lorentz force law 1)10. 12|l . but rather the dynamical 
variable p. This problem is simply addressed by replacing (|10.9() with 



Si^ j dX{i)iaj,ip+p{x-%lj^Qi}j)+qxA{x)) (10.14) 

so that the p and ij) equations become 

p^qF-i (10.15) 
ipia^ — pipjo- (10.16) 

The quantity p-x is a constant of the motion, and can be viewed as defining the 
mass. 

A more serious problem remains, however. If we form the spin bivector 
S = ijjia^il} we find that 

S^2pSx (10.17) 

so, if p and x are initially colincar, the spin bivector does not precess, even in 
the presence of an external S-field To solve this problem an extra term 

must be introduced into the action. The simplest modification is 

S2= JdX {ijia3^+p{x - i^-foij) + qxA{x) - ^Fi^ias^p) (10.18) 

which now yields the equations 

x^ipjo'4' (10.19) 

p = qF-x- —VF(x)-S (10.20) 
2m 

ijias = p^jo - -^Fil;. (10.21) 
Am 

The problem with this system of equations is that m has to be introduced 
explicitly, and there is nothing to identify this quantity with p-x. If we assume 
that p = mi we can recover the pair of equations 

S^—FxS (10.22) 



i ^ —F-v - -^WF(x)-S, (10.23) 

~" 2m^ 



which were studied in j73| . 

Whilst a satisfactory semiclassical mechanics for an electron still eludes us, it 
should be clear that the STA is a very useful tool in constructing and analysing 
candidate models. 

10.2 Grassmann Algebra 

Grassmann algebras play an essential role in many areas of modern quantum 
theory. However, nearly all calculations with Grassmann algebra can be per- 
formed more efficiently with geometric algebra. A set of quantities {Ci} form a 
Grassmann algebra if their product is totally antisymmetric 

= -OC (10.24) 
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Examples include fermion creation operators, the fermionic generators of a su- 
persymmetry algebra, and ghost fields in the path integral quantisation of non- 
abelian gauge theories. Any expression involving the Grassmann variables {Q} 
has a geometric algebra equivalent in which the {Q} are replaced by a frame of 
independent vectors {e^} and the Grassmann product is replaced by the outer 
(wedge) product |74l 175) . For example, we can make the replacement 

CiCj ^ SiAej. (10.25) 

This translation on its own clearly does not achieve a great deal, but the ge- 
ometric algebra form becomes more powerful when we consider the 'calculus' 
defined by Berezin [2]. This calculus is defined by the rules 



1. 

together with the 'graded Leibniz' rule', 



= (5,j, (10.26) 



Ql^=S,i. (10.27) 



^(/i/2) = |^/2 + (-l)[^^l/#, (10.28) 

where [/i] is the parity (even/odd) of /i. In geometric algebra, the operation of 
the Grassmann derivatives can be replaced by inner products of the reciprocal 
frame vectors 

A(^e^( (10.29) 

so that 

II - ^■e,=S]. (10.30) 

Some consequences of this translation procedure were discussed in j^, where 
it was shown that the geometric product made available by the geometric al- 
gebra formulation simplifies many computations. Applications discussed in |74| 
included 'Grauss' integrals, pseudoclassical mechanics, path integrals and Grass- 
mann- Fourier transforms. It was also shown that super-Lie algebras have a very 
simple representation within geometric algebra. There seems little doubt that 
the systematic replacement of Grassmann variables with geometric multivectors 
would considerably enhance our understanding of quantum field theory. 



11 Conclusions 

There is a growing realisation that geometric algebra provides a unified and pow- 
erful tool for the study of many areas of mathematics, physics and engineering. 
The underlying algebraic structure (Clifford algebra) appears in many key areas 
of physics and geometry [76) , and the geometric techniques are finding increasing 
application in areas as diverse as gravitation theory and robotics |77| [TSj . 
The only impediment to the wider adoption of geometric algebra appears to be 
physicists' understandable reluctance to adopt new techniques. We hope that 
the applications discussed in this paper make a convincing case for the use of ge- 
ometric algebra, and in particular the STA, in electron physics. Unfortunately, 
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in concentrating on a single area of physics, the unifying potential of geometric 
algebra does not necessarily come across. However, a brief look at other ap- 
plications should convince one of the wider utility of many of the techniques 
developed here. 

Further work in this field will centre on the multiparticle STA. At various 
points we have discussed using the multiparticle STA to analyse the non-locality 
revealed by EPR-type experiments. This is just one of many potential appli- 
cations of the approach outlined here. Others include following the streamlines 
for two particles through a scattering event, or using the 3-particle algebra to 
model pair creation. It will also be of considerable interest to develop simpli- 
fied techniques for handling more complicated many-body problems. Behind 
these goals lies the desire to construct an alternative to the current technique of 
fermionic field quantisation. The canonical anticommutation relations imposed 
there remain mysterious, despite forty years of discussion of the spin-statistics 
theorem. Elsewhere, there is still a clear need to develop the wavepacket ap- 
proach to tunnelling. This is true not only of fermions, but also of photons, on 
which most of the present experiments are performed. 

Looking further afield, the approach to the Dirac equation described in 
Section 0] extends simply to the case of a gravitational background ■ The 
wavepacket and multiparticle techniques developed here are essentially all that 
is required to address issues such as superradiance and pair creation by black 
holes. Closer to home, the STA is a powerful tool for classical relativistic physics. 
We dealt briefly with the construction of classical models for the electron in 
Section ^1 Elsewhere, similar techniques have been applied to the study of 
radiation reaction and the Lorentz-Dirac equation The range of applica- 
bility of geometric algebra is truly vast. We believe that all physicists should 
be exposed to its benefits and insights. 
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A The Spherical Monogenic Functions 

We begin by assuming that the spherical monogenic is an eigenstate of the a;AV 
and J3 operators, where all operators follow the conventions of Section We 
label this state as ip{l, so 

-xAVij{l, fi) = /?/-(/, fi), J^i^il, fi) = fiTpil, n). (A.l) 

The Ji operators satisfy 

JiJt = -[(icr,,)-(a;AV) - \iai][{ia.i)-{xhV) - \icTi] 

= 3/4-a;AV + (a;AVa;AV) (A.2) 

where the V indicates that the derivative acts on everything to its right. Since 
{xhVxhV)ij = xhV{xhVip) - xAVtjj (A. 3) 

we find that 

J^J^^P{l,^I) = [3/4+ 21 + 1^)^(1, fi) 

= (; + 1/2)(; + 3/2)^(/,m). (A.4) 

With the ladder operators J+ and J_ defined by 

J+ = Ji+jJ2 (A.5) 
J-=Ji-jJ2, (A.6) 



it is a simple matter to prove the following results: 

[ J+ , J_] = 2 J3 AJ, ^J-J+ + .h + ,h^ 

[J± , h] = TJ± J.Ji = J+J-- J3 + Ja^- 



(A.7) 



The raising operator J+ increases the eigenvalue of J3 by an integer. But, for 
fixed Z, /i must ultimately attain some maximum value. Denoting this value as 
IJL+, we must reach a state for which 

J+V;(Z,M+) = 0. (A.8) 

Acting on this state with Ji Ji and using one of the results in (|A.7|I we find that 

(/ + l/2)(Z + 3/2) = /i+(/i+ + l) (A.9) 
and, as I is positive and represents an upper bound, it follows that 

/i+=Z + l/2. (A.IO) 
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There must similarly be a lowest eigenvalue of J3 and a corresponding state 
with 

J-iJj{l,H-) = 0. (A.ll) 

In this case we find that 

(? + l/2)(? + 3/2) = /^_(/i_-l) (A.12) 

=^ fi-= -{1 + 1/2). (A.13) 

The spectrum of eigenvalues of J3 therefore ranges from (1 + 1/2) to —{1 + 1/2), 
a total of 2(/ + 1) states. Since the J3 eigenvalues are always of the form (integer 
+1/2), it is simpler to label the spherical monogenics with a pair of integers. 
We therefore write the spherical monogenics as ipj^, where 

-xAVi^r = ii^r i > (A.14) 

Jsi^r = {m+ Di^T -l-l<m<l. (A. 15) 

To find an explicit form for the we first construct the highest-m case. 
This satisfies 

J+^\ = (A.16) 
and it is not hard to see that this equation is solved by 

tA;' oc sin'6le"^'^^'^^ (A.17) 

Introducing a convenient factor, we write 

V^J = {21 + l)P/(cos0) e^'^^'^3_ (A.I8) 

Our convention for the associated Legendre polynomials follows Gradshteyn & 
Ryzhik j44|, so 

Pl^{x) = ^--^(1 - x2)"/2— — - 1)' (A.19) 



and we have the following recursion relations: 

(1 _ x2)^^n^ + mxPr{x) = -(1 - x^fl^Pr^\x) 



dx 



(A.20) 



The lowering operator J_ has the following effect on ^: 

J-iP = [-dei^ + cot6'90V'*f^3]e~'^*'^=' - ia2 ^{tp + crs^Paa). (A.22) 

The latter term just projects out the {l,icr3} terms and multiplies them by 
—ia2- This is the analog of the lowering matrix in the standard formalism. The 
derivatives acting on tpl form 

= {21 + l)[-dgPl{cos9) -lcotePl{cos9)]e^^ " l)'^^'^^ 

= {21 + l)2ZP/-i(cos6»)e(' " (A.23) 
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and, if we use the result that 

<J^ = <J2e't''''\ (A.24) 

we find that 

oc [2ZP/-i(cos^) - P/(cos^)iCT^]e(' " (A.25) 

Proceeding in this manner, we are led to the following formula for the spherical 
monogenics: 

^J" = [{l + rn + VjPricosO) - P^+^ (cos^)ic7^]e"*'^*^3 , (A.26) 

in which I is a positive integer or zero, m ranges from — (Z + 1) to I and the 
P™ arc taken to be zero if |m| >l. The positive- and negative-m states can be 
related using the result that 

P^^{x) = (-l)™|^p™(a;), (A.27) 

from which it can be shown that 

m-i<^2) = {-lr^^^^^i^'^^'^ . (A.28) 

The spherical monogenics presented here are unnormalised. Normalisation fact- 
ors are not hard to compute, and we find that 

d6 ^sinfiiVrV-r =47r^^77 (A.29) 

Jo Jo (' ~ '^j' 
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